Astronomy & Astrophysics manuscript no. overview 


© ESO 2009 


August 28, 2009 





C2D Spitzer-IRS spectra of disks around T Tauri stars 

IV. Crystalline silicates 

J. Olofsson', J.-C. Augereau^, E. F. van Dishoeck^'^, B. Merm^, F. Lahuis^'^, J. Kessler-Silacci^, C. P. Dullemond^, I. 
01iveira2-i2, G. A. Blake^ A. C. A. Boogert^ J. M. Brown^ N. J. Evans 11^, V. Geers"', C. Knez", J.-L. Monin\ and 

K. Pontoppidan^^ 

' Laboratoire d'Astrophysique de Grenoble, Universite Joseph Fourier, CNRS, UMR 5571, Grenoble, France 
e-mail: olof ssonOobs . u j f-grenoble . f r 

Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands 
Max Planck Institut fiir Extraterrestrische Physik, Giessenbachstrasse 1, 85748 Garching, Germany 
Herschel Science Centre, SRE-SDH, ESA RO. Box 78, 28691 Villanueva de la Caiiada, Madrid, Spain 
SRON Netherlands Institute for Space Research, P.O.Box 800, 9700 AV Groningen, The Netheriands 
The University of Texas at Austin, Department of Astronomy, 1 University Station C1400, Austin, Texas 787120259, USA 
Max Planck Institute for Astronomy, Konigstuhl 17, 69117 Heidelberg, Germany 

Division of Geological and Planetary Sciences 150-21, California Institute of Technology, Pasadena, CA 91 125, USA 
Infrared Processing and Analysis Center, MS 100-22, California Institute of Technology, Pasadena, CA 91 125, USA 
Department of Astronomy and Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON, M5S 3H4, Canada 
Department of Astronomy, University of Maryland, College Park, MD 20742, USA 

Division of Geological and Planetary Sciences 150-21, California Institute of Technology, Pasadena, CA 91 125, USA 



Received August 28, 2009; accepted 

ABSTRACT 

Aims. Dust grains in the planet forming regions around young stars are expected to be heavily processed due to coagulation, fragmen- 
tation and crystallization. This paper focuses on the crystalline silicate dust grains in protoplanetary disks, for a statistically significant 
number of TTauri stars (96). 

Methods. As part of the Cores to Disks (c2d) Legacy Program, we obtained more than a hundred Spitzer/IRS spectra of TTauri stars, 
over a spectral range of 5-35 fim where many silicate amorphous and crystalline solid-state features are present. At these wavelengths, 
observations probe the upper layers of accretion disks up to distances of a dozen AU from the central object. 

Results. More than 3/4 of our objects show at least one crystalline silicate emission feature that can be essentially attributed to 
Mg-rich silicates. Fe-rich crystalline silicates are largely absent in the c2d IRS spectra. The strength and detection frequency of the 
crystalline features seen ai A > 20yum correlate with each other, while they are largely uncorrelated with the observational properties 
of the amorphous silicate 10/jm feature. This supports the idea that the IRS spectra essentially probe two independent disk regions: 
a warm zone (<1 AU) emitting al A ~ lOyum and a much colder region emitting at A > 20 /jm (<10 AU). We identify a crystallinity 
paradox, as the long-wavelength {A > 20yum) crystalline silicate features are 3.5 times more frequently detected (~55% vs. ~15%) 
than the crystalline features arising from much warmer disk regions (A ~ 10/jm). This suggests that the disk has an inhomogeneous 
dust composition within ~10 AU. The analysis of the shape and strength of both the amorphous lOyum feature and the crystalline 
feature around 23 /jm provides evidence for the prevalence of /jm-sized (amorphous and crystalline) grains in upper layers of disks. 
Conclusions. The abundant crystalline silicates found far from their presumed formation regions suggests efficient outward radial 
transport mechanisms in the disks around TTauri disks. The presence of yum-sized grains in disk atmospheres, despite the short time- 
scales for settling to the midplane, suggests efficient (turbulent) vertical diffusion, likely accompanied by grain-grain fragmentation to 
balance the efficient growth expected. In this scenario, the depletion of submicron-sized grains in the upper layers of the disks points 
toward removal mechanisms such as stellar winds or radiation pressure. 

Key words. Stars: pre-main sequence - planetary systems: protoplanetary disks - circumstellar matter - Infrared: stars - Methods: 
statistical - Techniques: spectroscopic 



1. Introduction 

The silicate dust grains that are originally incorporated into stel- 
lar nebulae and eventually constitute planet forming disks are 
thought to be of interstellar medium (ISM) size and composi- 
tion, namely sub-micron in diameter and extremely amorphous 
(>99%) in structure (e.g. lGai]|ll998h . Solar System comets, on 
the other hand, show high crystallinity fractions. The silicate 
grains in Oort comet Hale-Bopp are for instance composed of 
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40 to 60% of Mg-rich cry stalline grains dWooden et al.lll999l 
120071 : ICrovisier et al.|[T997l) . while Jupiter Family comets have 
slightly lowe r crysta lline fractions (of about 25-35%, e.g. comet 
9P/Tempel 1, Ha rker e t al. 2007). Although the actual crystalline 
fraction may depend somewhat on th e methodology u sed to cal- 
culate the dust optical properties (e.g. lMin et al ] l2008h . the clear 
detection of crystalline silicates features in comet spectra indi- 
cates their mass fraction to be high compared to the ISM. The 
amorphous silicate grains that were incorporated into the Solar 
nebula must have been processed in the very early stages of 
planet formation before being incorporated into comets. Because 
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of the high temperatures required to thermally anneal amorphous 
silicates and modify their lattice structure, it is generally ac- 
cepted that the crystallization of amorphous silicates occured 
close to the young Sun. The crystalline silicates were then trans- 
ported outwards due to processes such as turbulence or winds, 
leading to possible radial variations that may explain some of the 
mineralogic differences observed between small Solar System 
bodies today. Other processes could also explain the presence 
of crystals in the outer regions, such as shock waves triggered 
by gravitational instabilit ies i n the solar nebula as de scribed in 
iDesch & Connollvl(l2002l) and lHarker & DeschI (|2002|) . Locally, 
such shocks could sufficiently heat up surrounding materials to 
anneal and crystallize amorphous grains. In all these scenarii, sil- 
icate crystalline grains can be considered as tracers of the history 
of the protoplanetary, circumsolar disk. 

The detection of silicates in the disks around young stars was 
notoriously difficult until the Infrared Space Observatory (ISO) 
and, later, the Spitzer Space Telescope were launched. Thanks to 
ISO, mid- to far-IR spectroscopy of disks around young interme- 
diate mass stars (Herbig Ae/Be stars, hereafter, HAeBe) could be 
performed. , Acke & van den Ancker (2004) showed that 52% of 
HAeBe stars exhibit emission in the 10 /im (amorphous) silicate 
feature that arises from the upper layers of the disks, while 23% 
of these stars show features at 1 1.3 jum, which is associated with 
forsterite, a Mg-rich crystalline silicate. Usin g the mid-IR ca- 
pabilit ies of MIDI on the VLT Interferometer. Ivan Boekel et alj 
( |2004|) spatially resolved the 10 //m emission zone for three 
HAeBe stars. They found a gradient of crystallinity, with a larger 
fraction of enstatite and forsterite crystalline grains located close 
to the star (1-2 AU) than further away in the disk (2-20 AU). 
This was among the first observations where the disk dust con- 
tent was proven to be radially inhomogeneous. 

Although TTauri sta rs were too faint f or observations with 
the ISO/SWS and LWS. lNatta etall ( l2000h employed the lower 
spectral resolution ISOPHOT instrument to detect the broad 
amorphous silicate emission feature at around 10 //m for 9 stars 
in the Chameleon I cloud; but the low resolution, the low 
signal-to-noise ratio and the limited spectral range did not al- 
low firm detection of any crystalline features in the disks around 
young solar mass stars. Pioneering ground-based observations 
of TTauri stars (hereafter TTs) nevertheless showed evidence 
for forsterite emissi on features i n the t errestrial spectral win- 
dow around 10 /vm. iHonda et al] ( l2003h . for example, clearly 
identified the presence of crystalline silicate emission features 
(forsterite at 10.1, 10.5 and 11.2yum, and possibly enstatite at 
10.9 jum) in the N-band spectrum {R ~ 250) of Hen 3-600A, 
and argued that the crystalline grains represent about 50% (30% 
enstatite, 20% forsterite) of the total mass of grains emitting 
at lOjum. Because of the difficulty of disentangling between 
grain growth, the emission from PAH (Polycyclic Aromatic 
Hydrocarbon) or from forsterite as the main cause of the 11.3 yum 
"feature" i n some disks (see e.g. ISitko et ani2000l and the dis- 
cussion in lPrzvgodda e t al. 2003b, an unambiguous assessment 
of the presence of crystalline silicates in large samples of TTauri 
disks has awaited the launch of Spitzer 

The high sensitivity and large spectral range (5-35 //m) 
of the IRS spectroscopic instrument onboard Spitzer permits 
routine mid-infrared spectroscopy of TTauri disks in nearby 
molecular clouds. Surveys of young s olar analogs i n diff erent 
star forming region s such as those by iFurlan et al.l (l2006l) and 
IWatson et al.l ( l2009h for the Taurus cloud, and of amorphous sil- 
icate features in other clouds by Kessler-Silacci et al. (2006), 
were performed. Detailed mineralogic fits analyzing the amor- 
phous and crystalline silicate fractions in TTauri disks were 



perfor med by iBouwman et al.l ( l2008i 7 TTs) and ISargent et al.l 
( l2009l 65 TTs), as well as for disks around much lower mass ob- 
je cts like the borderi ine brown dwarf (hereafter, BD) SST-Lup3- 
1 dMerm et al.ll2007[). and the BD 2M ASS J044427 13+25 12 164 
in the Taurus cloud ( Bouv et al.l2008l) . both of which show abun- 



dant crystalline silicate grains despite their low ( sub-)stellar tem- 
peratures. According to iBouwman et alj ( l2008h . the crystalline 
content of the inner regions of TTauri disks appear to be mostly 
dominated by enstatite silicates while the crystalline content of 
the outer region instead seems dominated by forsterite grains. 
Such data point toward an inhomogeneous dust content within 
the disks which could result from a radial dependence of crys- 
tallisation processes, or a difference in the initial conditions un- 
der which these crystals formed. The evolution of c rystalline 
silicates has also been studied by IWatson et al.l ( l2009l 84 TTs), 
who find that within the same environmental conditions, the sili- 
cate crystalline mass fraction varies greatly, from none to nearly 
100%. They find that 90% of the objects in their sample present 
crystalline features in the lOjum region (either the 9.2//m en- 
statite feature or the 11.3 jum forsterite feature). Regarding other 
crystalline features at longer wavelengths, they found that 50% 
of the objects surveyed had detectable crystalline features at 
~33 jum. Also, no strong correlations between crystalline indices 
and stellar parameters have been found, suggesting that some 
mechanisms (e.g. X-ray emission, giant planets formation or mi- 
gration) operate to erase the correlations that are expected from 
the standard models of crystallisation processes. 

Here we present a comprehensive study of crystalline sili- 
cates in more than a hundred disks around young stars based 
on 5-35 yum IRS spectroscopic observations obtained as part 
of the c 2d Spitzer Legacy Program From Molecular Clouds to 
Planets dEvans et al. I l2003h . This work is the continuation of a 
series of c2d/IRS papers, studying grains, PAHs and gas in in- 
ner disk regions around young stellar objects (class II). Paper 
I bv iKessler-Silacci et al.l (|2006|) focuses on amorphous silicate 
features as a proxy for grain growth in a subsample of the c2d 
IRS Class II objects. They f ind that the grain s ize is a function 
of spectral type. Paper II bv iGeers et al.l i2006l) shows evidence 
for low fractions of PAHs around TTs, indicating low gas phase 
PAHs abundance in disks as compared to the ISM. Paper III by 
iLahuis et al.l (l2007l) reports the first detections of [Ne II] and [Fe 
I] gas-phase lines in TTauri disks. 

In this paper, we present in Sec.|2]IRS observations of the en- 
tire c2d/IRS sample of class II objects, and we provide in Sec. [3] 
an overview of the detection statistics of solid-state emission 
features, especially those from crystalline silicates. In Sec.|4]we 
study the properties of crystalline silicate grains, including their 
size, and search for correlations between amorphous and crys- 
talline silicates and between the SED shape and grain size. We 
further investigate the difference in crystallinity between the in- 
ner and outer regions of disks. In Sec.|5]we discuss the implica- 
tions of our results for the dynamics of disks at planet-forming 
radii, and we summarize our results in Sec.|6l In a companion 
paper (Olofsson et al. 2009b, in prep.), we analyse in more de- 
tail the dust mineralogy based on compositional fitting of the 
complete spectra. 



2. Spitzer/IRS observations 

We present in this study the infrared spectra of disks around 108 
young stellar objects obtained as part of the c2d Spitzer Legacy 
program. The spectra were obtained using the IRS instrument 
onboard the Spitzer Space Telescope. It was used to expand the 
early spectroscopic studies of HAeBe stars with ISO, and of a 
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few TTs observed from the ground in the A ~ 10 /^m atmospheric 
window. Thanks to the high sensitivity of Spitzer, a large sam- 
ple of disks around low-mass and solar-mass stars could be effi- 
ciently observed. 

The IRS spectrograph is composed of four different modules 
that enable a wavelength coverage from 5 to 37 fim. Two mod- 
ules, Short-Low and Long-Low (SL and LL hereafter), provide a 
spectral resolution between 60-127 over ranges of 5.2 - 14.5 jum 
and 14.0-38.0/im for SL and LL, respectively. The remaining 
two modules, Short-High and Long-High (SH and LH hereafter), 
have a spectral resolution of about 600. They cover more limited 
spectral ranges of 9.9- 19.6;um and 18.7 - 37.2yL(m for SH and 
LH, respectively. 

2.1. Stellar sample 

In this paper, we focus on Class II objects - young stars no 
longer embedded in their protostellar envelope and surrounded 
by a circumstellar disk, where the dust may already be signif- 
icantly processed. Our selection contains objects classified as 
Class II according to the literature, but also 10 objects with 
no known spectral type. In this sample, 30 objects are sources 
newly discovered with the instrument IRAC onboard Spitzer as 
part of the mapping phase of the c2d project. All such objects 
are named with the p refix SSTc2d, and are classified as Class 
II objects according to lEvans et al. I (120091) except for one source 
with no known classification (SSTc2d J182909.8+03446). The 
object SSTc2d J 16222 1.0-230403 is characterized as Class I, but 
shows amorphous 10 /im emission without additional absorp- 
tion features. Two other SSTc2d objects have no classification. 
However, for the 10 objects with no classification in the litera- 
ture, they were retained in the full list as they clearly show amor- 
phous silicate features in emission. 

The final sample contains 108 stars including 60 TTs, 9 
HAeBe and 1 BD. The 38 remaining objects have no spec- 
tral classification in the literature, or are the new SSTc2d stars 
discovered with Spitzer. This sample is distributed toward sev- 
eral clouds: Perseus (16 objects), Taurus (9), Chamaeleon (23), 
Ophiuchus (25), Lupus (16), Serpens (15), and includes 4 iso- 
lated stars: 3 HAeBes (BFOri, HD 132947 and HD 135344) and 
one TT star (IRAS 08267-3336). All six clouds are young star 
forming regions (1-5 Myr), located at distances within 140 and 
260 pc from the Sun, with star densities (number of young stellar 
objects per square parsec) lying between 3 and 13. Some clouds 
are therefore quite extended , like Taurus or Lu pus. The main 
cloud properties are given in lEvans et al.l ( |2009[) . The complete 
list of IRS targets analysed in this paper, their classifications and 
spectral types, as well as details about the observations, are given 
in Tablegl 

Out of the 108 objects, four display spectra unusual 
for Class II objects. These objects are DLCha, HD 132947, 
SSTc2d J 182849.4+00604 and SSTc2d J182859.5+03003. The 
first three sources all display spectra reminiscent of those ob- 
served toward Ox ygen-rich Asymptotic Giant Branch stars (see 
ISloan et al1l2003l for example spectra) while the fourth object 
displays a spectrum similar to a C-rich AGB star. Both SSTc2d 
objects located in the Serpens cloud are objects newly discov- 
ered using Spitzer/IRAC surveys. DLCha is a known variable 
star in Chamaeleon and HD 132947 is an isolated Herbig star, 
classified in the literature as a Pre-Main Sequence star (see e.g. 
[Valenti et al. 2003). We do not question here the classification of 
these four stars but as they display unusual spectra (see Fig.l24]i 
we decided to remove them from the rest of the study. 



2.2. Data reduction 



The spectra presented in this paper were extracted from S13 
pre-reduced (BCD ) data using the c2d legacy team pipeline 
dLahuis et al ] l2006h which makes available two different extrac- 
tion methods: full aperture extraction and optimal Point Spread 
Function (PSF) extraction. The full aperture extraction method 
is applied to both low and high spectral resolution modules. For 
the low resolution modules, the c2d pipeline implements an ex- 
traction with a fixed-width aperture over the whole order The 
source position in the slit is determined and the aperture is then 
centered on the source. The width is such that 99% of the flux 
of a point-source falls within the window. For the high resolu- 
tion modules, the full slit width is used for the extraction. One 
disadvantage of this method is that the reduced spectra present 
a number of spikes (mainly due to bad pixels response, or "hot" 
pixels on the array, see below) that cannot be easily removed by 
hands-off pipelines. 

For the PSF extraction, the observed signal is assumed to be 
that of a point source or slightly resolved source plus an uniform 
zero level. This zero level, in most cases, represents the local 
extended emission close to the source, and the typical PSF pro- 
file is calculated from high signal-to-noise ratio calibrator data. 
As this zero-level may also contain residuals from, for example, 
raw pipeline dark current, we remove this contribution. PSF fit- 
ting is less sensitive to bad data samples and unidentified bad 
pixels than the aperture extraction method. The Spitzer Science 
Center provides masks with known bad pixels, but they are not 
all identified. Our pipeline detects such remaining pixels. The 
LH array is particularly affected by this issue, but other modules 
require attention. With the PSF extraction method, bad pixels are 
simply eliminated, whereas their values are simply interpolated 
for full aperture extraction method. The PSF extraction method, 
in addition, provides an estimate of the data zero level or sky 
contribution to the observed spectrum. An important disadvan- 
tage of the PSF fitting method, though, is that for some modules 
(mainly for the 2"'' orders), the PSF is subpixel in size and the 
extraction method can become unstable. Thus, for some sources, 
the 2"'* order of the SL module and the 2"'' order of the LL mod- 
ule are unavailable, leading to incomplete spectra. 

In this paper, we adopt the PSF extraction method in or- 
der to obtain spectra with as few spikes as possible. Where the 
PSF extraction method is not available or turns out to be unsta- 
ble, we use the full aperture extraction method to build spec- 
tra. We then opt for a median filtering scheme to remove the 
remaining spikes: for the SL and SH modules we use a median 
smoothing over three channels while we use a median smooth- 
ing over five channels for the LL and LH modules. Some prob- 
lems were encountered with the data reduction for two objects in 
our sample: EC69 and ISO-Cha237. For the latter source there 
are two possible explanations: either the LLl module spectrum 
poorly matches the LL2 module spectrum, or the LL2 module 
pipeline reduction failed and its slope is off. The second expla- 
nation seems the most likely, given the shape of the red side of 
the lOyum feature. For EC69, flux values for some orders were 
null. We decided to keep these two objects in our sample, first 
because ISO-Cha237 shows a 10 /^m feature that can be ana- 
lyzed. For EC69, the spectrum appears featureless, and as it is 
classified as a TTauri star, we choose not to bias our results by 
removed "problematic" spectra. 
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Table 1. Overview of the modules used in all the different cases 
encountered to obtain the final spectra (see Sec. |2.3| for details). 
The exact modules used for each individual star are given in 
Tables 



Table 2. Summary of the different crystalline silicate features 
examined in this study. 



Case 


SLl 


SL2 


SH 


LLl 


LL2 


LH 


Number of stars 









V 






V 


9 


1 






V 


V 




V 


2 


2 


V 


V 




V 


V 




51 


3 


V 


V 


V 






V 


18 


4 


V 


V 


V 


V 




V 


22 


5 


V 


V 


V 


V 


V 




6 



Note: SLl: 7.4-14.5 //m , SL2: 5.2- 
LL2: 14.0-21.3/im 



7.7;um,LLl: 19.5 - 38.0A'm, 



2.3. Module merging and offset correction 

The c2d pipeline corrects for telescope pointing errors that 
may lead to important flux offsets between the different mod- 
ules. Some small offsets can nevertheless remain between mod- 
ules even after pointing correction, and an additionnal post- 
processing offset correction is thus applied. This correction de- 
pends on the modules available for every spectrum. Table[T]gives 
a summary of all the possible configurations we encountered. 
A tick-mark means that data provided by this module are 
used to build the final spectrum. Priority was given to SL and 
LL, first because the spectral resolution is sufficient for detect- 
ing crystalline features, and second, because the high resolution 
modules provide noisier spectra than those from the low reso- 
lution modules. For each star, a final spectrum is then obtained 
by merging all the chosen modules and correcting manually for 
small remaining offsets. The spectra themselves are presented in 
Fig.ffl Fig.lHand Figs.[l5H23] 

2.4. Estimating the uncertainties 

The uncertainties propagated through the c2d pipeline are domi- 
nated by the flux calibration error estimates. As our study is ori- 
ented toward the detection of features in emission from silicates, 
we are only concerned with relative uncertainties which are eval- 
uated by dividing the original spectrum by a s moothed version of 
itself obtained using a Savitzky-Golay filter (ISavitzky & Golayl 
Il964l) . This provides an estimate of the signal-to-noise ratio 
(SNR) on the spectrum, and is used to assess the presence of 
features as follows: The Savitzky-Golay filtering process calcu- 
lates a local 3"'-order polynomial for every point of the spec- 
trum, based on the four left and right neighbors. The main ad- 
vantage of this filter is that it maintains the shape and contrast of 
the features, while "pointy" features tend to be rounded out with 
classical average smoothing, thereby increasing the uncertainties 
at the top of every feature and downgrading the actual quality of 
the data. The Savitzky-Golay filter has only been used to calcu- 
late the uncertainties displayed as grey envelopes in Fig.IU Fig. El 
and Figs. [T5l - l23l not affecting the spectra themself presented on 
these figures and used for the scientific analysis. 



3. Overview of observed solid-state features 

3.1. Observed crystalline silicates 

The 104 spectra analyzed in this paper are displayed in 
Figures [TS] to l23l for wavelengths between 5 and 35 fim and in 



Crystal 


Peak 


Measured 


Measured AA 




A [jjm] 


peak A [jjm] 


Eq. Width [/jm] 




9.2"''^* 


9.24±0.08 


0.49±0.11 


Enstatite 
(MgSiOs) 


21.6' ' 
23.0<'-2> 
24.5' • ' 


21.47±0.16 

C23 

Cz3 


0.45±0.25 

C23 

C/3 




28.2(i-2) 


C28 


C28 




11.3*^' 


11.26±0.12 


0.60±0.33 


Forsterite 

(Mg2Si04) 


16.2<2) 
23.8<2' 
27.6<2) 


16.08±0.16 

C23 

C28 


0.38±0.16 

C23 

C28 




33.6<2) 


33.1+0.63 


1.19±0.47 


Diopside 
(CaMgSiaOe) 


25.0'2' 


25.14±0.18 


0.58±0.22 



References: " iKoike et al1 ( l200(]h . '^ iMolster et afl ( |2002|) 



units of Jy (Fy). They are sorted by cloud and increasing Right 
Ascension. The spectra of the four isolated stars in our sample 
are displayed in Fig. [23] 

Broad emission features at 10 and 2Qfim from amorphous 
silicates can be identified in many spectra, as well as nar- 
rower emission features. The latter are more easily identified 
beyond 20yum, and are due to crystalline silicates. Figure[T] 
illustrates the diversity of the spectra seen in our sample. 
ISO-ChaII54 (bottom spectrum) shows many crystalUne fea- 
tures, with a weak lOjum band. In contrast, AS 205 (top spec- 
trum) is almost "pristine", with strong amorphous bands and 
no evidence of crystalline feature^ The third spectrum, for 
SSTc2d J033036.0+303024, has both amorphous and crystalhne 
features. The crystalline features in disks can, in most cases, be 
attributed to Mg-rich silicate minerals (with a possible contri- 
bution from Ca-rich silicates); in agreement with some model 
predictions (e.g.iGail 2004 and ref. therein), a nd as observed in 
solar system comets (e.g. lWooden et al]|2007h . 

3.1.1. Mg-rich crystalline silicates: Forsterite & Enstatite 

Table|2] presents the different crystalline silicate features we 
searched for in every spectrum - in addition to the amorphous 
silicate feature at lOyum. In Table|2] "Peak A" is taken from 
the literature, "Measured /}" and "A/l" values are the peak po- 
sitions and equivalent widths measured by our observations. 
"C23" and "C28" refer to silicate features blended in complexes. 
Uncertainties are standard deviations. Here we consider fea- 
tures arising from Mg-rich silicates, specifically enstatite (inosil- 
icate belonging to the pyroxene group, with chemical formula 
MgSiOj) and forsterite, the nesosilicate Mg-rich end-member of 
the olivine group (Mg2Si04). Blowups of representative spectral 
features are shown in Fig.|2] Using IRS spectra, we cannot dis- 
entangle differences between crystalline structures, for instance, 
clino-e nstatite or ortho-enstatite. Still, according to Koike et id] 
(I2000h . enstatite features are most likely attributed to ortho- 
enstatite. 

The example spectrum for ISO-Chall 54, displayed in Fig.[T] 
as well as the blowups in the right panel of Fig.|2] shows that 
the 23.0, 24.5/^m enstatite and the 23.8yum forsterite features 



' Extensive, narrow features not discussed in this paper, but seen 
between 10 and 2 /im in Fig.[T] are attributed to water emission lines 



l lSalvk et alj|2008l) 
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0.01 I , , , , I lJ . I , I d , , , I , L, I I . I I , , LJ , I , , , 

5 10 15 20 25 30 35 

A [fJ-m] 



Fig. 1. Example Spitzer/IRS spectra of 3 TTauri stars (top spectrum: SL+SH+LH modules, bottom spectrum; SL+LL). The grey 
envelopes correspond to 3-cr uncertainties as estimated in Sec. 12. 41 The vertical lines show the peak positions of enstatite (dashed 
lines) and forsterite (dot-dashed lines) crystalline features that were searched for in every spectrum of our 108 star sample (see also 
Table|2]i. The 12.5jum feature (triple-dot dashed line) is attributed to silica, and the two yellow boxes correspond to the C23 and 
C28 crystalline complexes discussed in Sec. 13.1.11 



can be blended into one single complex. Similarly, the 27.6yum 
forsterite feature can also be blended with the 28.2yum enstatite 
feature. In the following, we will treat these features as two 
complexes, at reference positions of 23 and 28 /^m, indepen- 
dent of the actual crystals responsible for the observed emission. 
We shall note these two complexes, C23 and C28, which are 
by far the most frequently detected crystalline features in disks 
(Sec. |3.2.TT i. Fig.|2]displays eight example spectra of the C23 and 
C28 complexes. Clearly, the shape and peak positions may vary 
from one object to another We will show in Sec. l4.1l that some 
of these differences may relate to the mean size of the emitting 
grains, although the composition (more or less enstatite com- 
pared to forsterite), grain shape, and strengh of the continuum 
might also contribute to the diversity of the C23 and C28 fea- 
tures observed. 

Beside the features listed in Table|2] other enstatite bands at 
10.6, 11.7 and 19.6/^m, as well as forsterite features at 10.1 and 
19.7 fim were also searched for with little success in a subsample 
of 47 high SNR spectra. The detection rate of these crystalline 
features is very low and not significant. For the two enstatite 
and forsterite features at ~19.6 and 19.7yum, the main difficulty 
lies in the added noise where SH and LH modules merge, which 
prevents any firm detection of these two features. 

3.1.2. Fe-rich components: Troilite & Fayalite 

The formation of Mg-rich crystalline silicates by direct con- 
densation from a gas of solar composition (as opposed to for- 
mation by thermal annealing) is expected to additionally yield 
iron-bearing condensates, in particular FeS (or troilite, see 
[Wooden et al. 2007). This condensate is of great interest and, 
while a detailed study of the ISO spectra of t wo HAeBe stars (in- 
cluding HD 163296, Fig.[T8]l carried out bv lKeller et all (l2002b 



indicates FeS in these systems, it remains difficult to detect, in 
general, because the main feature (between 17-26 /^m), centered 
at 23.5 //m, is broad, and falls in regions where the IRS spectra 
are often noisy (due to module merging around 20 /im). Further, 
the disk spectra are often dominated by the presence of the 1 8 //m 
amorphous silicate feature as well as by frequent C23 emission 
features. Therefore, studying the presence or absence of FeS 
emission in the c2d sample is beyond the scope of this paper 

Another Fe-rich component that can potentially be detected 
in IRS spectra is fayalite, a crystal belo nging to the olivine 
class (Fe2Si04). Theoretical opacities from lFabian et al.l (1200 ll) 
show strong emission features around 30 ;um. For small grains, 
three narrow features are expected at ~27, 29.3 and 31.6yL(m. 
We find there are 4 objects where fayalite may tentatively be 
present, all four in the Chameleon cloud (Fig.fTTTi: TWCha, 
VZCha, WXCha, and Sz62. Interestingly, the Sz62 spectrum 
also shows other features that can be attributed to fayalite, at 
roughly 11 //m and 19/im. All these objects notably display a 
feature centered at around 30.5 //m, which is at shorter wave- 
length comp ared to the feature f ro m pure fayaUte w hich lies at 
31.6yum (e.g. lFabian et al.ll200ll or lKoike et al.ll2003l) . This may 
reflect a small departure from pure Fe-rich crystalline olivine 
(e.g. ~10-15% Mg fraction) as illus trated by the three top spec- 
tra in Fig. 1 from lKoike et al.l (l2003l) where it is seen that the fea- 
ture shifts to smaller wavelengths with increasing Mg-content. 
We note that there are other spectra that show some similar be- 
haviours in the 30 fim region, but they are of lower signal-to- 
noise ratio and contain possible data reduction artifacts (see e.g. 
RXJ0432.8-H1735 spectrum, Fig.[T6ll. These factors make it dif- 
ficult to conclude whether fayalite is present or not. Because of 
the small number of positive detections and because of the am- 
biguity between possible noise and real features, we will not in- 
vestigate the presence of this component further. We do note, 
however, that the compositional fits to seven Spitzer/IRS spectra 
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Fig. 2. Le/f panel: Examples of the 9.2 (enstatite), 11.3 (forsterite) and 12.5 yum (silica) features superposed on the amorphous lOjum 
feature. Right panel: Examples of the C23 and C28 crystalline silicate complexes. The grey areas correspond to 3-cr uncertainties. 
Dashed lines correspond to enstatite emission features, dot-dashed lines correspond to forsterite emission features. The sources in 
the left panel are a) WaOph6, b) LkHQ'325, c) TWCha, d) SSTc2dJ161 159.8-382338, e) SSTc2d J033036.0-H303024, f) Sz52, 
g) XX Cha and h) ROX42C. 



of TTs by iBouwman et alj ( l2008h show no evidence of Fe-rich 
materials. 

However, the previous conclusion means that we do not ob- 
serve pure fayalite in our spectra. Considering the following for- 
mula for crystalline olivine (Mg vFei_^)2Si04 (with < jc < 1), 
features in the 20-30//m spectral range become weaker (espe- 
cially around 23 //m) for Mg fractions below 0.6. This means 
that we are actually probing crystalline olivine with at least more 
than 60% of Mg compared to Fe. Still, as crystalline pyroxene is 
also contributing to C23 and C28 complexes, further interpreta- 
tions remain difficult. 



Finally, silica (SiOi) has also been identified in our spec- 
tra, with a feature arising at 12.5yum. Attribution of this sin- 
gle fe ature to amorphous o r crystalline silica is not straightfor- 
ward. ISargent et al.l (l2009l) showed that amorphous silica tends 
to produce a feature around 12.3-12.4yum, while various crys- 
talline silica polymorphs produce a feature centered at or slightly 
to longer wavelengths of 12.5 //m. Due to the presence of the 
broad amorphous lOjum feature, disentangling the different con- 
tributions is beyond the scope of our statistical study of the 
most prominent crystalline silicates. Example silica features are 
shown on Fig.[T]and the left panel of Fig.|2l 



3.1 .3. Other silicates: Diopside & Silica 

Diopside, a calcium magnesium silicate, and a member of the 
pyroxene group (CaMg(Si03)2), can be searched for at ~25 jum. 
Because of the presence of the C23 complex, however, the diop- 
side feature can sometimes be blended with the complex. Even 
worse, the shoulder of the complex can be mis-interpreted as a 
25 jum feature. We therefore obtain only a lower limit on the fre- 
quency at which this feature is present. The low detection rate 
of an additional (but weaker) 20.6 [im diopside feature does not 
strengthen the confidence level of the detection of the diopside 
in TTauri star disks. 



3.2. Silicate feature statistics for TTauri stars 

3.2.1 . The fraction of disks showing silicate(s) feature(s) 

We developed a routine that measures the characteristics of both 
crystalline and amorphous (lOjum) features. For each feature, 
it assumes a local continuum which is built using a two point 
Lagrange polynomial joining the feet of the feature, computed 
as follows: 
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where ^ is the final continuum, Ai and A2 the blue and red feet 
of the feature, respectively, and Fy, and Fy,, the observed fluxes 
at A] and A2, respectively. This routine returns the band flux, its 
peak position, the mean wavelength and the peak flux. It also 
returns the l-cr uncertainty on the band flux. The results given 
by this procedure are presented in Table|5] in the form of the 
SNR on the measured feature fluxes using the above procedure. 
Results for only 101 objects out of total 104 are shown since 
3 of them (2.8% of total sample) do not display any crystalline 
or amorphous silicate features. Those objects are: EC 69, EC 92 
and Sz 84. 

Thus, the overwhelming majority (~97%) of the disks do 
show at least one feature that can be attributed to silicates. Since 
10//m-sized or larger grains are essentially featureless in the 
IRS spectral range (see, for example, the theoretical opacities 
in Fig.|9l), our results qualitatively indicate that inner disk at- 
mospheres are populated by grains a few micrometers or less 
in size, independent of the age of our stars, and of cloud mem- 
bership. Some differences between clouds are nevertheless ad- 
dressed in Sec. 13.2.21 and a discussion concerning the size of 
the amorphous and crystalline silicate grains can be found in 
Sec.O 

In the following, we remove the 8 HAeBe stars (plus 
HD 132947 that we already removed) in our survey from the sta- 
tistical analysis so that the sample is mainly composed of objects 
classified as TTs. The statistics in the paper are therefore made 
on the remaining 96 objects. 

3.2.2. The prevalence of crystalline features 

Feature-by-feature statistics are displayed in Figure[3] which 
shows detection fractions of crystalline features ranging between 
~10% and ~60%. A visual inspection of the 28-29 yum and 33- 
35 jum features for a subset (40 TTs and 7 HAeBes) of th e targe t 
sample analyzed in this paper led iKessler-Silacci et all (|2006|) 
to conclude that ~50% of the spectra show crystalline silicate 
features. Considering now the full 96 star sample and only se- 
cure detections with SNR > 20 for the C28 and 33.6//m com- 
plexes, we obtain a fraction of ~65% (63 out of 96) for ob- 
jects displaying at least one of the two features. Including all 
the other crystalline features (still with SNR > 20), this fraction 
rises to ~78% (75 out of 96 stars). We have checked that these 
results are not corrupted with possible PAH emission features, 
and that the lL3um can be associated with forsterite. Indeed, 
iGeers et al] (l2006h identified only three objects in our sample 
that have potential PAHs emission at 11.3 yu m: LkHQ-330 , TCha 
and SR21. All the other objects cited in jGeersetal] (|2006|) 
are HAeBe stars and , therefore, no t considered in our sample. 
With additional data, iGeersI ( l2007h could confirm the presence 
of PAHs for LkHa 330, T Cha and SR 21, and so these three ob- 
jects have then been removed from our statistics for the anal- 
ysis of the 1 1. 3 yum forsterit e feature. One other object, EC 82 
is identified bv iGeersI (|2007|) to have a 6.2 yum PAH feature but 
with no 1 1.3 fim counterpart, thereby avoiding any possible con- 
tamination in this particular case. Note that we have been as 
conservative as possible in identifying features, especially for 
the crystalline features around lOjum. Some spectra display a 
shoulder at around 11.3jum that could possibly be interpreted 
as a forsterite feature (see e.g. Sz52 spectrum, Fig.|2] label/). 
We choose not to attribute such kind of shoulder to forsterite. 
Therefore, we conclude from our analysis that about 3/4 of the 
protoplanetary disks in our sample of young solar-analogs show 
at least one crystalline feature in their 5-35 yum spectrum, mean- 




9.2 10.0 11.3 12.5 16.2 21.6 C23 25.0 C28 33.6 



Fig. 3. Detection statistics for crystalline silicate features (yellow 
and red bars) for the lO.O/zm amorphous silicate feature (blue) 
and for silica at 12.5 yum (green). Darker color bars are positive 
detections (SNR > 20), and lighter color bars stand for tenta- 
tive detections (SNR < 20). The uncertainties are those due to 
Poisson statistics. 



ing that dust crystallization is not a marginal process but is in- 
stead widespread at this stage of disk evolution. 

It is interestin g to compare these results to those of 
IWatson et al.l (l2009l) who analyzed 84 TTs in the Taurus-Auriga 
association, also using Spitzer/IRS. Developing indices mea- 
suring the departure from an ISM-like pristine lOyum feature, 
they found that 90% of the Taurus-Auriga disks show either the 
9.2 jum enstatite feature or the 11.3jum forsterite feature, while 
we ob tain about 20% for both features (Fig.|3]l. In lWatson et al.l 
(I2009h . any departure from the lOyum pristine reference spec- 
trum is considered as being evidence for crystalline features, 
while we argue that this frequent departure that we also observe 
has rather more to do with a grain size effect (see Sec. l4.1l and 
Sec. |5.3| for further discussions). There is therefore a difference 
in the interpretati on of the mea s ureme nts that l eads to inconsis- 
t ent co nclusions. I Watson et al.l (l2009l) used the i D'Alessio et all 
(2006) model to investigate the impact of dust sedimentation on 
the observed dispersion in equivalent width for the lOyum amor- 
phous feature. They could reproduce a broad range of observed 
equivalent widths by depleting the dust fraction (that consists 
of grains < 0.25 yum in their model) in the upper layers of the 
disks, thereby mimicking dus t settling. In a more recent paper, 
iDullemond & DominikI ( |2008|) find that sedimentation has, nev- 
ertheless, a very limited effect on the shape of the lOyum feature 
as long as a continuous grain size distribution is considered (as 
opposed to fixed grain sizes). For longer wave length crystall ine 
features, comparison with the work of iWatson et al.l (|2009|) is 
possible for the 33.6yum forsterite feature, but not for the fea- 
tures around 23 and 28 yum as they did not study these two com- 
plexes. They find that 50% of the Taurus-Auriga objects show 
the 33.6yum crystalline feature, compared to 28% for our sam- 
ple, but this detection statistic is consistent with the statistics for 
the C23 (55%) and C28 (54%) crystalhne complexes. 

Figure|4|shows feature-by-feature statistics sorted by clouds. 
The largest occurence of the C23 complex is found in the 
Ophiuchus cloud with about 80% of the objects showing this 
feature. Disks in Serpens have a prevalence of crystalline fea- 
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tures significantly smaller than in the Perseus, Chamaeleon, 
Ophiuchus and Lupus clouds which show comparable fractions 
of amorphous and crystalline features. The Taurus cloud shows a 
noticeably low frequency of the lOjjm amorphous feature com- 
pared to the other clouds and does not appear to be balanced 
by a remarkably high fraction of crystalline grains when com- 
pared to the other regions. Our sample contains only eight ob- 
jects hiJheJ[aurusclou^ however, which introduces a strong 
bias. IWatson et al.l ( l2009l) found that the amorphous 10 //m fea- 
ture is present for a large majority of the objects in their survey 
of 84 TTs. Therefore, discussing the differences is not statis- 
tically significant for our Taurus sample. One can notice that, 
except for Taurus, the five remaining clouds show very simi- 
lar behavior. This points toward a rather homogeneous distri- 
bution of amorphous and crystalline silicate dust grains in disks 
throughout clouds. Detection statistics in the Serpens cloud are 
overall smaller, but the trends are the same as that in the other 
clouds. Our current result on the Serpens cloud will be revisited 
in a forthcoming paper by Oliveira et al. (in prep) who analyze 
a much larger sample of more than a hundred Serpens objects 
observed with Spitzer/IRS. 

In the following, we will treat all the objects as a single sam- 
ple. The first reason for this is that cloud-to-cloud differences 
suffer from low statistics and merging the entire sample will in- 
crease the significance of our inferences (96 objects in total). 
The second reason is that the trends of the different clouds seem 
to be similar (except for our Taurus cloud smaple which has a 
low number of stars). Merging the sample will strenghten our 
conclusions thanks to the large number of Class II objects. 



4. Properties of crystalline silicates in the disks 
around young stars 

One may expect that the relatively large spectral range of the 
IRS instrument (5-35 jum) implies that the detected short- and 
long-wavelength emission fea tures could arise from phys ically 
distinct regions i n the disk. iKessler-Silacci et al] ( l2006h and 
iPinte et al.r ( l2008h . among others, have shown that for a typi- 
cal TTs, the 10 //m emission comes from a region about 10 times 
smaller than the 20 jum emitting zone. We examine in this section 
the correlations between the detected features. In particular, we 
inspect the relation between the properties of the 10 //m silicate 
amorphous feature (its energy, the size of its carriers, its appari- 
tion frequency) and those of the crystalline silicate features at 
A > 20jum. 



4. 1. Growth of warm amorphous silicates 

The amorphou s lOyum feature is present at a rate of about 70% 
in our sample. iBouwman et al" (12001 1) . Ivan Boekel et al.1 ( l2003h 
and lKessler-Silacci et al. 1 (120061) have shown that detailed studies 
of this feature reveal much information on dust characteristics. 
Especially noteworthy is that the observed shape versus strength 
relation for the 10/im feature can be reproduced by varying the 
characteristic grain size over the range of 0.1-3 //m. Here we 
analyse the lO/^m feature by adopting the same continuum nor- 
malization and the same computational method for the shape 



(flux ratio S m/^Sg g) and strength (Sp^J^^) of the lOyum feature 
as in lKessler-Silacci et al. I (120061) . The different values 5 11.3, 5 9.8 
and 5 Peak are obtained by normalizing the observed flux Fy as 
follows: S'v = 1 + {Fy - Fy^c)/{Fv^c}, where Fy^^ is the estimated 
linear continuum between the red and blue feet of the feature and 
(Fy^c) its mean value between the two feet (all fluxes expressed 
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Fig. 4. Same as Fig. [3] but now sorted by clouds. Uncertainties 
are again Poisson noise, derived from the number of objects per 
cloud: 15 for Perseus, 8 for Taurus, 20 for Chamaeleon, 16 for 
Lupus, 24 for Ophiuchus and 12 for Serpens. 



in Jy). We choose a wavelength range of +0.1 fxm around 9.8yL/m 
and 11.3 /im to calculate 5 1 1 .3 and 5 9.8 . 

Figure|5] (left panel) shows the known correlation between 
the flux ratio 5 n. 3/59.8 (lO/^m feature shape) and S'peak (lO/^m 
feature strength). To quantify this correlation we use a Kendall t 
test. The Kendall r rank correlation coefficient measures the de- 
gree of correspondance between two datasets. If the agreement 
between the two sets of values is perfect then r = 1, if the dis- 
agreement is perfect then t = -I. The Kendall t procedure also 
returns a probabilit y P, computed usi ng the complementary error 
function of \t\ (see lPress et al.|[T992l for further informations on 
the Kendall t and P calculations). The P probability is the two- 
sided significance; the smaller the P value, the more significant 
the correlation. 

For the shape versus strength of the amorphous 10/im fea- 
ture, we find an anti-correlation (r - -0.42) with a signif- 
icance probability P = 3.58 x 10"^. This correlation can be 
interpreted as larger grains producing a flatter lOfim amor- 
phous feat ure, while smaller gra ins produce narrower fea- 
tures (e.g. van Boekel et ap |2003|). A thorough anal y sis by 
iKessler-Silacci et alj (|2006|) and iKessler-Silacci et all (120071) 
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of the shape-strength trend showed that more massive stars 
(HAeBe) tend to have sharper 10 //m features than do TTs and 
BDs, which cluster in a region more consistent with micron- 
diameter grains as the typical emitting grain size. The compar- 
ison of our observations to synthetic 10 jum amorphous features 
calculated for different grain sizes (see Sec. l4.6l for details on the 
computation method of grain opacities) confirms this trend as 
the bulk of our sample clusters in a region consistent with grains 
that have grown to sizes larger than about 2 fim (left panel of 
Fig.lH). 

4.2. Growth of cold crystalline silicates 

Although a similar analysis of the IQ/im amorphous emission 
features is rendered more difficult because of the broadness 
of the feature, the narrower longer-wavelength crystalline sili- 
cate features, especially the C23 complex, show a grain size- 
dependent behaviour that can be tested with our disk sample (see 
the theoretical opacities in Fig.|9]l. Similar trends could not be 
exploited for the C28 complex as its shape is not much sensitive 
to grain size (over the 0.1-3 yum range). For the C23 complex, we 
build a feature shape index, the S24/S23 flux ratio, computed us- 
ing the same approach as for the amorphous IQfim feature. The 
5 23 normalized flux is obtained by averaging the normalized (lo- 
cal) continuum-subtracted flux between 23.5 and 23.8 ;um, and 
between 24.1 and 24.4 //m for 524. Figure|5] (right panel) shows 
the S24/S23 ratio as a function of the strength of the C23 com- 
plex (C23 5 Peak)- The Kendall test gives t - -0.54 with a sig- 
nificance probability below 10"^*^, indicative of a clear anticor- 
relation between these two quantities. 

To better evaluate the relationship between grain size and 
the C23 shape, we calculate theoretical S24/S23 ratios for syn- 
thetic crystalline silicates opacity curves obtaine d using th e 
Distribution of Hollow Spheres (DHS) method (Min et al .120051) . 
We consider grains between 0.1 and 2.7 /im in radius, and two 
different comp ositions: 100% forsterite (Mg2Si04, using optical 
constants from lServoin & Pirioull 19731) and a mixture composed 
of 50% forste rite p lus 50% enstatite (MgSiOj, optical constants 
taken from Jaeger e t al. 1998). Grains larger than ~ 3.0 /urn, for 
either enstatite or forsterite, do not give striking results, as the 
feature disappears into the continuum. Also, a pure enstatite 
composition turns out to be incompatible with the observed 
shape for the C23 complex as it shows two unblended features at 
23.8 and 24.5 i-im. This indicates that we are not observing pure 
enstatite grains in our sample, a result w hich seems to first order 
consistent with iBouwman et al. 1 (I2008h who found that the very 
inner (~ 1 AU) warm dust population is dominated by enstatite 
while the dust population at large radii is dominated by forsterite. 
The overall shape-strength trend for the C23 crystalline complex 
is better reproduced for a forsterite/enstatite mixture, with char- 
acteristic grain sizes larger than 1.0//m in order to to match the 
bulk of the measured values. This suggests that not only have the 
warm amorphous silicate grains grown in our sample of TTs (left 
panel of Fig.jSj, but also that the colder crystalline silicate grains 
have coagulated into /.<m-sized particles and remain suspended 
in the upper layers of the disk atmospheres. 

The offset in the right panel of Fig.|5]between the data points 
distribution and the curves obtained using theoretical opacities 
may come from various effects. First, we consider pure crys- 
talline dust grains, which may be a too restrictive assumption 
to reproduce the observations. It allows us to see the impact of 
grain size, however, and to sh ow that it matches reasonably well 
the observed trend. Recently, iMin et al" I (120081) studied the im- 
pact of inclusions of small crystals into larger amorphous grains. 
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Fig. 6. Left panel: Correlations between the slope of spectra 
{F30/F13) as a function of the strength {S peat) of the amor- 
phous 10/^m feature. The two outliers (blue open circles) are 
cold disks: LkHQ'330 and CoKuTau/4. Right panel: Spectral 
slope {F-iolFi-},) versus strength (S'peak) for the C23jum crys- 
talline complex. 

but this is too much of a refinement for our purposes. Second, 
we only consider pu re Mg-rich crystals (no iron). According to 
IChihara et alj (|2002|) . the 24.5 jum feature of crystalline pyrox- 
enes shifts towards slightly la rger wavelengths \y ith increasing 
Fe-content. Also, according to lKoike et al.l (l2003b . the peak po- 
sition of the 23.8/im crystalline olivine feature does not shift 
unless the Fe-content is larger than 20%. As the 524 index is 
computed integrating between 24.1 and 24.4//m, raising the 
Fe-content will produce a smaller 5 24 and therefore a smaller 
S24/S23, while the strength of the complex will slightly decrease 
as the 24.5 jum feature shifts to larger wavelengths. To conclude, 
even if there is a small shift when comparing the observations 
to theoretical points, the trend is well reproduced by varying the 
grain size. 

Interestingly, we find no correlation between the grain size 
proxy for the amorphous 10 pm feature and the grain size 
proxy for the C23 complex. Considering the sets of values 
(5ii.3/59.8)/5^^'^i, as a function of {S24/S23)/Sll^^, we indeed 
obtain a Kendall's t of 0.07, with a significance probability P 
of 0.44. This means that even if grain growth took place in both 
components, the growth of the amorphous and crystalline grains, 
respectively, does not seem to be correlated. 

4.3. Relationship between the 10pm and C23 features with 
disk properties 

Young flared disks are expected to show rising spectra in the 
mid-IR, as opposed to more evolved, settled systems having 
self-shadowed disks with much flatter mid-IR spectra (e.g. 
iDominik et al.ll2003h . Here we measure the slope of the spectra 
by the ratio of the F30 and F13 indexes, where is the mean 
flux value (in Jy) in the range 13//m±0.5 pm, and between 29 pm 
and 31 pm for Fjo. The left panel of Figure |6] shows the corre- 
lation between the strength of the amorphous lOyum feature and 
the slope of the spectra (F30/F13 ratio). Taking out the two out- 
liers (blue open c ircles), wh ich are known cold disks (LkHo' 330 
and CoKu Tau/4, iBrown et al. 2007), we find a correlation co- 
efficient r = 0.29 with a significance probability P smaller than 
5.2 X 10""* suggestive of a correlation. A similar trend is found by 
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Shape (5'24/5'23) versus strength (S'peak) for the C23 ;um crystalline complex. Colored points and curves are derived from theoretical 
opacities for different mixtures. For every curve on left panel, the open circles represent different grain sizes, from bottom to top, 
0.1, 1.0, 1.25, 1.5, 2.0, 2.7, 3.25, 4.3, 5.2, and 6.25 /^m. For every curve on right panel, the grain sizes are, from bottom to top, 0.1, 
1.0, 1.25, 1.5, 2.0 and 2.7 yum. 



iBouwnian et al.' ('2008") within their FEPS sample of seven stars, 
and by IWatson et al.. (2009.) with the Taurus-Auriga sample of 



4.4. Relationship between ttie 10 /im feature and crystallinity 



84 TTs. This can be interpreted as smaller grains (large S^pJ^^ 
values) being probed in flared disks (high values of /^3o/^^i3), 
while flattened disk (low F^q/Fu values) atmospheres are char- 
acterized by larger grains. 

If dust sedimentation is the main process responsible for disk 
flattening, one would expect the upper layers of the disks to be 
depleted in large grains, hence the mid-infrared spectroscopic 
signature to be that of very small grains (peaked 10//m features). 
The observed trend (Fig.|6]l very clearly refutes such a possi- 
ble trend in the relation between the SED shape and that of the 
10/^m feature. The effect of dust sedimentation on the shape of 
the lOyum amorphous silicate feature was recently discussed by 
[OuUemond & Dominik (2008) based on models. They showed 
that the sedimentation of the larger grains toward the disk mid- 
plane is generally not enough to transform a flat lOfim feature 
into a peaked feature. Our observations tend to support these re- 
sults as sedimentation, if assumed to be revealed by disk flatten- 
ing, is generally associated with flat 10//m features. This trend is 
in fact valid for the majority of the TTs in our sample as the ob- 
jects cluster in the region with flat SEDs and flat silicate features 
in Figure |6] 

For the C23 crystalline complex (right panel of Fig.|6]l, we 
find a rather weak anticorrelation, with r = -0.21 and a sig- 
nificance probability P <0.01. The small crystalline grains with 
5 Peak larger than about 1.2 are only found in flat or slowly ris- 
ing disks (F30/F13 smaller than about 2), while the larger grains 
can be found in all kinds of disks. Interestingly, the very flared 
disks which have large crystalline grains emitting at 23 fim also 
show small amorphous grains in their warm component (that is, 
the region probed at 10/im). Seven out of the eight objects (RY 
Lup, RXJ1615. 3-3255, Haro 1-1, SSTc2d J162148.5-234027, 
SSTc2d J 1 62245 .4-243 1 24, SSTc2d J 1 627 1 5 . 1 -245 1 39, Haro 1 - 
16) whith Fio/Fii > 3 on the right panel of Fig.|6]indeed have 
their counterpart in the left panel with larger than about 2 
(the last object being DoAr 25). 



We find in Sec. 14. 21 that the growth of the warm amorphous 
silicates (probed by the lOjum feature) is statistically discon- 
nected from the growth of the cold crystalline silicates (probed 
by the C23 complex). We further investigate the relationship 
between the amorphous and crystalline silicates features by 
searching for a correlation between the energy contained in the 
amorphous lOyum feature and the energy contained in the C23, 
C28 and the 33.6 ;um features. To get rid of possible distance 
or brightness effects we normalized the line fluxes as follows: 
Line flux [unitless] = Line flux[W.m"^]/(/l(F^_c)), where (F.i.c) 
is the mean value of the local feature continuum and A is the 
measured central wavelength of the considered feature. 

The unitless line fluxes are displayed in the left panel of 
Figure|2l which shows that when both the C23 and the 10 //m 
feature are present with SNR > 20, their emission energies are 
not correlated (t - 0.06, significance probability P - 0.55), al- 
though the plot would suggest that objects with low lOfim line 
flux usually also have a low C23 line flux. The middle panel of 
Figure|2]shows a similar result when both the C28 and the lOjum 
feature are present with SNR > 20: r = 0.023 (no correlation be- 
tween these features) with a significance probability P - 0.83, 
reflecting the dispersion. Finally, the right panel of Fig. [T] shows 
a similar trend for the 33.6 jum and the amorphous lO/^m feature: 
T - 0.22 and a significance probability P - 0.15. Therefore, the 
flux radiated by the 10 //m amorphous feature is largely unre- 
lated to the energy contained in the crystalline silicates features 
that appear at wavelengths longer than 20 /im. 

Similarly, the shape of the amorphous 10/im feature (that 
is, the 5 11. 3/59.8 ratio, which is a proxy for grain size. Sec. 14. It 
seems uncorrected with the emergence of the 23 //m complex, 
the 28 //m complex, or the 33.6//m forsterite feature. The right 
panel of Fig.|2] for example, shows the number of positive de- 
tections of the C23, C28 and 33.6//m features per 5 n j/S'y.g bin, 
where the errors bars correspond to Poisson noise (square root 
of the number of detections). The fraction of disks with C23 and 
C28 features is found to be independent of the mean grain size 
probed by the 10 //m amorphous feature, further suggesting that 
the lOfim feature and the C23, C28 features essentially probe 
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Fig. 7. Left panel: From left to right: correlations between the unitless line fluxes of the \Q^m feature and the C23, C28 complexes 
plus 33.6 fim forsterite feature. Right panel: From left to right: percentages of objects showing the C23/C28 complexes and 33.6 //m 
forsterite feature as a function of the shape of the IQpm feature. Larger, amorphous grains are located at the top of the plot. 



disconnected populations. For the 33.6yum feature, the fraction 
of objects per bin seems less independent because of the smaller 
numbers of detections (23 objects with both the 10 and 33.6^(m 
features). 

Overall, these results suggest that the crystalline features ap- 
pearing at wavelengths longer than 20 //m and the 10 //m amor- 
phous feature arise from unrelated populations. On the other 
hand, Figure[8]shows the correlation between the normalized line 
fluxes of the C23 and C28 complexes where both are present 
with SNR > 20, and also their correlation with the 33.6/im fea- 
ture. The Kendall t tests and associated significance probabil- 
ities confirm the expected result that the crystalline features at 
wavelengths longer than 20 pn\ have emission fluxes that corre- 
late to each other. 

4.5. Feature correlation coefficients 

The one-to-one correlation probabilities between the detection 
frequencies of all the features examined in this study with SNR 
> 20 are displayed in Table[3] The diagonal shows the total num- 
ber of positive detections for any given feature (these numbers 
being close to percentages of detections as we have 96 objects). 
The lower part of the Table displays the Kendall t values as well 
as the number of spectra with common detections of the consid- 
ered features between parenthesis. The Kendall significance P 
probabilities are displayed in the upper part. 

At this point, it is interesting to understand the meaning of 
the inferred Kendall t values in this Table. Some correlations 
can possibly be found according to the Kendall t values but 
only because both features are absent from spectra. Consider 
as an example the 11.3/im forsterite feature and the C23 com- 
plex, for which we obtain t = 0.244 and P - 0.001, for 12 
detections of the 11.3 //m feature and 56 detections of the C23 
complex, indicating a weak correlation. Out of the 96 stars, 1 1 
spectra show both features, one spectrum shows the 11.3 fim fea- 
ture only, and 45 objects show the C23 complex only. The other 
39 (96 - 45 - 1 - 1 1) objects show neither the 1 1.3 /im feature 
nor the C23 complex, and these objects account for concordant 
pairs. This explains why we find a weak correlation (non-zero t 
value), as the number of concordant pairs (50) is of the order of 
the number of discordant pairs (46) for this example. In case we 
would eliminate from the statistics objects that do not show any 
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Fig. 8. Top panel: Normalized fine fluxes for the 33.6 fim feature 
and the C23 complex (Kendal r - 0.55 and significance proba- 
biUty P = 2.41 x lO""*). Middle panel: Same but with the C28 
complex (t - 0.56, P = 3.66x 10"^). Bottom panel: Normalized 
line fluxes between the C23 and C28 complexes (t - 0.54, 
P = 1.19 X 10-^). 



of these two features, we would obtain an anti-correlation with 
T = -0.25 and P - 0.004. Therefore, in order to have compa- 
rable numbers for all the features, irrespective of the number of 
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Table 3. Correlation coefficients between the features detected with SNR> 20. Entries below the diagonal give the Kendall t values 
and the number of spectra with common detections of considered features between parenthesis. The diagonal displays the total 
number of detections for a given feature, while the entries above the diagnoal report the Kendall P significance probabilities. 
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detections, we decide to consider the entire sample for the cal- 
culation of the correlation coefficients, even when spectra do not 
show any of ffie considered features. 

The detection frequencies of the lOyum amorphous feature 
and the C23 complex are found to be essentially uncorrelated 
(t - 0.051 with P = 0.47). Similar results are obtained when 
considering the lO/im feature and the C28 complex (r = 0.013 
with P = 0.86), and also the lO/zm and 33.6yum features 
(r = 0.055 with P = 0.44). For the C23, C28 complexes and 
the 33.6/^m feature, we have, respectively, 43, 40 and 22 com- 
mon detections with the amorphous 10 /im feature (which is ob- 
served in 69 spectra). The C23 and C28 complexes, as well as 
the 33.6 fim forsterite feature, do therefore appear in spectra in- 
dependently from the presence of the lOyum feature. On the other 
hand, the C23 complex appears to be coincident with the C28 
complex in our spectra (r - 0.619 with P < 10"^^), which con- 
curs with our results based on line fluxes (Sec. |4.4t . The 33.6 fim 
feature also appears to be coincident, to a slightly lesser extent 
with the C23 complex (r = 0.288 with P < 10"^**), and the C28 
features (t - 0.233 with P - 3x 10^^). The number of common 
detections also confirm ffiese trends. 

Table [3] further shows that the forsterite 11.3yum feature is 
anti-correlated with the amorphous lOjum feature (r = -0.44, 
only 3 objects show both the lOyum and the 11.3/im features). 
The anti-correlation between the forsterite 11.3jum feature and 
the amorphous lOjum feature could be explained by a contrast 
effect. The amorphous feature is indeed rather stronger than the 
11.3 fim feature which could therefore potentially be present but 
hidden in the amorphous feature (see discussion in Sec. |4.6| for a 
more detailed analysis of this effect). 

Interestingly, the crystalline features in the 10 fim range (9.2 
and 11.3 fim) seem to be correlated with both the C23 and C28 
complexes: if crystalline features are present in the 10 fim spec- 
tral range then the presence of the two complexes is highly prob- 
able. Still the reverse consideration is not true, if one of the two 
complexes (C23 or C28) is present, it is not guaranteed that we 
will find ffie 9.2 or 11.3 fim features. This reflects that the cor- 
relation is dominated by objects where none of the considered 
features is present. Furthermore, as we find no strong correla- 
tion for the line fluxes regarding these features, this would mean 
that there is a general increase of crystallinity for several objects. 
Also, the two features at 9.2 and 1 1.3 fim do not seem to be cor- 
related with the 33.6 fim forsterite feature. The declining quality 
of the data (lower SNR) at the very end of the IRS spectra may 
explain this trend in part, as some (lower contrast) 33.6yum fea- 
tures may not have been detected. Another possible explanation 



would be that the 33.6 fim forsterite feature may probe an even 
colder disk region than the population probed by C23 and C28 
complexes. 

To summarize, our results show that the C23, C28 and 
33.6 fim crystalline features probably arise from the same pop- 
ulation of grains. Also, the detections of the 9.2 and 11.3jum 
crystalline features are correlated with the detections of C23 and 
C28 (but the presence of one of the longer wavelengths com- 
plexes does not involve the detection of 9.2 or 1 1 .3 fim features), 
which can be interpreted as a general increase of crystallinity in 
the inner and outer regions of disks. However, this crystalline 
population seems to be disjoint from the population of amor- 
phous sihcate grains that produce the 10 fim feature. 

4.6. Degree of crystallinity of the warm versus cold disk 
regions 

Both Figures [3] and |4] show larger occurrences for the C23 
and C28 emission features than for smaller wavelength crys- 
talline silicate features, especially the 9.2 (enstatite) and 1 1.3 fim 
(forsterite) features. As silicates are expected to thermally anneal 
close to the star, this result may reflect a counterintuitive positive 
gradient of the crystalline silicate fraction toward colder temper- 
atures. The presence of the broad, and generally strong, amor- 
phous lOyum feature may lower the contrast of (or even hide) the 
shortest wavelength crystalline features, however, preventing di- 
rect conclusions. We examine in this section this apparent crys- 
tallinity paradox by comparing our Spitzer/IRS observations to 
synthetic spectra. 

To better gauge the actual differences between the degrees of 
crystallinity of the warm {A ~ 10 fim) versus cold (A > 20jum) 
disk regions, we generate synthetic spectra representative of the 
observations, and from which we extract line fluxes as was per- 
formed for the IRS spectra in Sec. 13. 2.11 We consider amorphous 
species that include silicates of o livine stoichiometry (g lassy 
MgFeSi04, optical constants from iDorschner et al.l 1 1 9951) . and 
silicates of pyroxene stoichiometry (glassy MgFeSiiOe, opti- 
cal constants from Dorschner et al. 1995). For the crystalline 
species, we use enstatite (]V[gSi03) optical constants from 
iJaeger et"an (|1998[) and forst erite (Mg2Si04) optical constants 
from Servoin & Pirioul(ll973h . Theoretical mass opacities ka are 
computed using Mie theory (valid for hard spheres) for the amor- 
phous species an d DHS theory (Distribution of Hollow Spheres, 
iMin et all l2005h for the crystalline silicates. Twelve example 
opacity curves for the four compositions and three typical (in 
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Fig. 9. Example mass opacities (in arbitrary units) used to model the IRS observations. From bottom to top: amorphous olivine, 
amorphous pyroxene, crystalline forsterite and crystalline enstatite. The left panel shows opacities for a grain size a - Q.I fim, 
a = 1.5 fim for the middle panel and a - 6.0//m for the right panel. Some of the features discussed in Secs. l3.1.Tl and [3T31 can be 
recognized. Yellow boxes are places of both C23 and C28 crystalline complexes. 




Fig. 10. Synthetic spectra, in arbitrary units, for mass fractions of crystalline silicates from 0% to 50% ( steps are: 0, 5, 10, 15, 
20, 25, 30, 40, 50%, see Sec. 14. 61 for more details). The impact of the grain size on the detectability of the 11.3yum, C23 and C28 
features is considered for 0.1 /^m, 1.5 /^m and 6.0 fim grains. 



terms of spectroscopic signature) grain sizes (0. 1 jjm, 1 .5 fjm and 
6.0 fim) are displayed in Fig.|9] 

The synthetic spectra of Fig. [10] have been generated by 
adding to a representative continuum (Cv), the mass opacities 
of amorphous and crystalline silicates, multiplied by relative 
masses and by a blackbody at reference temperatures (350 K for 
the warm region, and 100 K for cold region). The amorphous 
content is a 50:50 mixture of olivine and pyroxene and the crys- 
talline content is a 50:50 mixture of enstatite plus forsterite. The 
adopted lO^m representative continuum is given by the median 
values (mean slopes and mean offsets: Cy = axA+b) of all the lo- 
cal continua linearly (CJ, - a' xA + b') estimated in Sec. |3.2.1| for 



the observed amorphous 10/im features: Cy - 0.01 x A + 0.097, 
with Cy in Jy and A in fim. Similarly, a representative C23/C28 
continuum has been derived for cold synthetic spectra: Cy = 
-0.001 xA + 0.43 (Cy in Jy and A in fim). 

To obtain as representative synthetic spectra as possible, we 
compute the median observed continuum over feature ratio (i.e., 
the continuum flux divided by the height of the feature) for the 
10/im feature and the C23 complex. For the amorphous lOyum 
feature, we obtain a ratio of 2 and we therefore normalize the 
opacities, before adding them to the continuum (Cy), so that the 
synthetic spectra show similar ratios. For the C23, we measure a 
median continuum over feature ratio of 20±10, and we normal- 
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Fig. 11. Signal-to-noise ratio of the 9.2jum, 11.3//m, C23 and C28 integrated line fluxes measured as in Sec. l3.2.ll but on the 
synthetic spectra displayed in Fig.[TO](see text for details about the models). The light gray area represents the conservative detection 
threshold used in this paper (SNR=20), while the dark gray area shows the SNR=3 limit. Grey diagonal stripes are the loci where 
theoretical opacities are not consistent with typical observations. 



ize the opacities in order to get a ratio between 10 and 30 in the 
synthetic spectra. This exercise could not be reproduced for the 
large 6.0 fim grains because the shapes of their C23/C28 features 
are not representative of the observed features (especially for the 
forsterite grains, see right panel of Fig.|9]l, suggesting such large 
grains are not the main carriers of these features. As a final step, 
the synthetic spectra are degraded by adding random noise, with 
a maximum amplitude of 1 .8 x lO""' Jy for the warm spectra and 
7.2 X 10"^ for the cold spectra, representative of the quadratic 
sum of the measured uncertainties on the Spitzer data over in 
the lOyum region and in the C23/C28 spectral range. The syn- 
thetic observations displayed in Fig. [10] are then processed as in 
Sec. |3.2.1| to derive the band fluxes of the crystalline features. 

The theoretical SNRs for the 9.2, 11.3//m, C23 and C28 
crystalline silicate features, assuming crystalline mass fractions 
between 0% and 50% and characteristic grain sizes, are dis- 
played in Figure[TT] It is seen that the SNR reached in the sim- 
ulations compare well with the measured SNR values, espe- 
cially for micron-sized grains (Table|5l). The diagonal-striped ar- 
eas correspond to crystalline mass fractions where the contin- 
uum over C23/C28 feature ratios could hardly match the median 
observed values, or which totally fail to reproduce them. This 
means that we are not able to discuss the detectability of features 
for crystalline mass fractions smaller than 5 and 15% for grains 
sizes of 0. 1 and 1 .5 /urn, respectively. In the particular case of 
the 6.0jum grains, as explained above, we cannot compute SNR 
for C23 and C28 complexes as opacities do not match typical 
observations. 

We see in Fig.[TT] that the SNR of the crystalline features 
drops with increasing grain size (especially the 9.2 jum, 1 1.3 yum 
features) and with decreasing crystalline mass fraction, as ex- 
pected. A crystalline fraction of ~ 5% or larger seems to be 
enough such that all four of the studied crystalline features are 
observable if the grains are sub-micron in size. For the 1.5 fim- 
sized grains, crystalline mass fractions of about 15% are re- 
quired, (1) to obtain spectra consistent with observations, and 
(2) to get the all the features studied detected with SNR>20. 
For higher crystalline mass fractions, they all present very simi- 
lar behaviors. For 6.0 //m-sized grains, it becomes more difficult 
to trace the detectability of the C23 and C28 features as opaci- 
ties become much less consistent with typical observations. For 
9.2yum and 11.3yum, crystalline fractions larger than 20% and 



15%, respectively, are required in order to detect these features 
with SNR>20. 

Since the observed features in our sample are mainly pro- 
duced by jum-sized grains (Sec. 14. 11 ), which in particular pro- 
duce a shoulder at the location of the 11.3 fim crystalline feature, 
we conclude that the physical processes that preferentially place 
jum-sized grains in the disk atmospheres do not contribute much 
to the measured differences between the number of detections of 
the 11.3 jim feature compared to the number of detections of the 
C23 complex (Figs.[3]and|4|. 

To summarize, we find an apparent crystallinity paradox, 
namely a noticeable difference in detection statistics for the 
11.3/im and C23 features, with more than 3 times more de- 
tections for C23 compared to the number of detections of the 
1 1.3 ^m. forsterite feature. This seems counterintuitive as shorter 
wavelengths are expected to probe warmer disk regions where 
grains can more efficiently crystallize. We investigated the pos- 
sibility that crystalline features could be hidden by amorphous 
features. The effect seems insufficient to explain the observa- 
tions, but we cannot yet firmly predict the real impact of this 
effect until the Spitzer observations are confronted by a more 
detailed compositional analysis. This is deferred to a future pa- 
per (Olofsson et al. 2009b, in prep.). 

5. Discussion 

5.1. On the lack of Fe-rich silicates 

Because of the large number of objects analysed in this pa- 
per, our results raise challenging questions on the dynamics 
and (chemical) evolution of planet forming disks regions that 
are in principle not statistically biased. For instance, the lack 
of Fe- rich silicates in our 96 star sampl^l reinforces the re- 
sult by iBouwman et al.l (|2008|) obtained on their 7-star sample. 
Silicates crystals studied in this analysis may contain iron, as 
seen in Sec |3.1.2| but magnesium is the most abundant element. 
Several scenarios can explain this apparent lack of iron. First, 
iDavoisne et al.l (l2006l) found that a reduction reaction during 
the thermal annealing of ferro-magnesian amorphous silicates at 
temperatures below lOOOK can produce pure forsterite crystals 

- except 5 marginal detections discussed in Sec. l3.1.2l 
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Fig. 12. Correlation between shape and strength for the amor- 
phous lOyL/m feature compared to models assuming a power law 
grain size distributions with indexes between -4 and -3. For each 
colored curve, from right to left, open circles correspond to mim- 
imum grain sizes equal to 0.1, 0.5, 0.75, 1.0, 1.5, 2.0 and 
4.0 fim. For blue points only the last four mininum grain sizes 
are represented, for red points the 0. 1 //m point is not shown, as 
their ^i'^'f" is larger than 3.5. 



Fig. 13. Correlation between shape and strength for the amor- 
phous lOyum feature compared to models assuming a variation 
in crystallinity fraction, for three different grain sizes. For each 
grain size, increasing crystallinity goes from bottom left to top 
right, from 0% to 50% crystalHnity (0, 5, 10, 15, 20, 25, 30, 40, 
50%) 



plus spheroidal metallic particles. The iron is therefore locked 
inside a metallic nanophase that cannot be observed. Secondly, 
according to Nuth & Johnson (2006), simple thermal anneal- 
ing in the warm inner regions of the disk cannot produce Fe- 
rich silicate crystals. As the annealing times cale for iron sili- 
cates is larger than for magnesium silicates (iHallenbeck et al.l 
I2000I) . the lifetime of such Fe-rich silicates against the evapora- 
tion timescale is expected to be too short for these grains to crys- 
tallize and then be transported outward in the disk. Still, thermal 
annealing induced by shocks may be able to produce both Fe- 
rich and Mg-rich crystals at the same time. 

5.2. On the need for turbulent diffusion and grain-grain 
fragmentation 

The new grain size proxy introduced here for cold crystalline 
silicates (via the C23 complex) indicates that the upper layers 
of the disks are preferentially populated by yum-sized grains. 
This result also holds for the warm disk region using the clas- 
sical size proxy for the carriers of the 10/im amorphous fea- 
ture (Fig.|5j. Adopting a purely sequential evolutionary approach 
for the solid constituents in disks, where interstellar-like grains 
coagulate to form micron-sized dust aggregates, then pebbles 
and planetesimals, a straighforward conclusion from these re- 
sults would be that grain agglomeration occured in the first few 
AU and that we are witnessing an intermediate, regular step 
of the growth of solid particles along their way to form plan- 
ets. Nevertheless, this picture is most certainly an oversimplified 
view of the processes at work in planet forming regions around 
young stars as it does not fit with basic dynamical considera- 
tions. Indeed, the short settling time-scale of jum-sized grains in 
a laminar disk, and t heir expected fast growth which will accel- 
erate sedimentation ( iDullemond & Dominikl 120051 : iLaibe et al] 
l2008h . both argue in favor of some processes that sustain the 
disk atmosphere with grains that have mid-infared spectroscopic 
signatures (grains smaller than about 10 //m) over a few Myr. 



Therefore, our results support the replenishment of the disk 
upper regions by the vertical transpor t of d ust particles (e.g. 
turbulent diffusion, Fromang & Nelsonl l2009l) . likely combined 
with grain-grain fragmentation to balance the expected effi- 
cient growth of solid particulates (that accelerate s sedimentation, 
IDullemond & Dominikll2005l: iLaibe et al.ll2008l) . In this picture, 
the grains observed at mid-IR wavelengths would then result 
from the destruction of much larger solid particles (an assump- 
tion consistent with the presence of millimeter-siz ed grains in 
most TTauri disks, e.g. review bv lNatta et al.ll2007l) . rather than 
from the direct growth from interstellar-like grains. 

5.3. On fhe depletion of submicron-sized grains 

A consequence of the above conclusion is that one would have to 
explain the general apparent depletion of submicron-sized grains 
in the disk zones probed by mid-lR spectroscopy. The fact that 
we see a flat, boxy lOyum feature profile for most TTs, indeed 
indicates that there cannot be many submicron-sized grains, be- 
cause their emission would overwhelm the few large grains. 
Grain-grain collisions are nevertheless expected to produce frag- 
ments with a range of sizes that are likely to extend to the submi- 
cron domain (i.e. the size of the monomers that form /^m-sized 
aggregates). This effect might be enhanced by the large velocity 
fluctuations (a fraction of the local sound speed) found at high 
altitudes ab ove the disk midplane in d isk simulations with MHD 
turbulence dFromang & Nelsonl2009l) . Qualitatively, submicron- 
sized grains, lifted in the disk atmospheres by turbulent diffu- 
sion or produced locally, are therefore expected to be abundant, 
which contradicts the IRS observations. 

5.3.1 . Impact of the size distribution 

The above conclusion is based on our single, characteristic grain 
size approach (Secs. |4n and |4.2t and one may question whether 
this assumption affects our interpretation of the observations. To 
answer that question, we calculate mean cross-sections as fol- 
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lows: 

(Tabs = I •^fl^Gabs(a) dn(a)/ I dn(a), 

where a is the grain radius, Qabs is the dimensionless Mie emis- 
sion/absorption coefficient for a 50:50 mixture of amorphous py- 
roxene and olivine, and dn{a) oc a''da is the differential grain 
size distribution. We arbitrarily choose Omax = 100/^m (large 
enough such that it does not affect the calculations for MRN- 
like size distributions), and we vary the minimum grain size 
between 0.1 and A-.Qfim. Using a classical MRN {p = -3.5) 
grain size distribution, we find that the diversity of 5 1 1 .3 /S' g.g and 
'^Peak" ^^li^ss is well reproduced by adjusting the minimum grain 
size flmin, and that the bulk of the observations is reproduced for 
Gmin larger than 2 fim (Fig.fTSli. a result consistent with our single 
grain size analysis. Nevertheless, another way to avoid having 
peaked lOyum amorphous features typical of pristine grains is to 
use a flatter slope for the size distribution. For example, p values 
larger than -3 reduce the contribution of submicron-sized grains 
to the mean absorption cross-section (see Fig.[T2]with p = -3), 
creating an amorphous 10 jum feature that is large and boxy de- 
spite the presence of some submicron-sized grains. Therefore, 
the majority of the points is located at the left side of Fig.[T2l 
which either probes yum-sized grains as the small end of the size 
distribution, or much flatter grain size distributions compared to 
the MRN, meaning that the amount of small grains is at least 
severely diminished in the regions probed by the Spitzer/IRS. 

5.3.2. Impact of the crystallinity 

Seve ral studies (e.g. B ouwman et ani200lLlHonda et al.ll2006l or 
Ivan Boekel et al. 2005) mention the possible impact of the crys- 
talUnity on the shape versus strength correlation for the lOyum 
feature. According to these studies, the prese nce of cry s talline 
grains can mimic the flattening of the feature. iMin et al.l ( l2008h 
recently investigated in detail this effect by considering inho- 
mogeneous aggregates. Their Fig. 2 shows that Spl^^ remains 
approximately constant for crystalline fractions up to 20% (the 
maximum fraction considered in their study) while 5 1 1.3 /S' 9.8 
rises, suggesting a shift along the y-axis of Fig. [12] due to crys- 
tallinity rather than along the observed trend. 

In order to further investigate the impact of crystallinity 
on this correlation, we use the synthetic spectra described in 
Sec. 14. 61 The only difference is that we do not normalize opaci- 
ties anymore to obtain a constant ratio between feature strength 
and continuum. Normalizing opacities would lead to a constant 
'^Peak"' which is not what we want to characterize here. We there- 
fore normalize the opacities, for each grain size, so that only the 
strongest feature reaches a ratio of 2 between the continuum and 
the feature strength. Other features are normalized with the same 
factor and obviously present a higher ratio between continuum 
and feature strength. This normalization is done in order to re- 
produce spectra similar to observations, and still conserve the di- 
versity of feature strengths. Subsequently, values for 5 9.8, 5 11.3 
and 'S'pII^j™ are derived in the same way as in Sec |4.1| Fig.[T3] 
shows the impact of crystallinity on the shape versus strength 
for the lOyum feature. We consider three grain sizes (0.1, 1.5 and 
6.0yum) and the same composition as in Sec. 14. 61 50:50 mixture 
of olivine and pyroxene for the amorphous content and 50:50 of 
enstatite and forsterite for the crystalline population. The crys- 
talUnity fraction varies from 0% to 50% (0, 5, 10, 15, 20, 25, 30, 
40 and 50%), and increasing crystallinity goes from bottom left 
to top right on the plot. We note that for 6.0 jum-sized grains the 



effect is negligible, while it is more dramatic for 0. 1 jum-sized 
grains. The spread induced by an increase of crystallinity is al- 
most orthogonal to the direction of the correla tion. This result 
was also found bv iKessler-Silacci et alj (|2006|) (Fig. lOe). The 
direction of this deviation can simply be understood when con- 
sidering opacities for amorphous and crystalline grains. For the 
same mass of dust, crystalline grain emission will be stronger 
compared to amorphous grains. Therefore adding some crys- 
talline grains to a purely amorphous dust content will give a 
lOyum feature stronger and broader (because of crystalline fea- 
tures at, e.g., 11.3 yum in the case of forsterite). This will thus pro- 
duce a higher value of 5 11.3/59.8 and at the same time a stronger 

■^Peak ■ 

Crystallinity therefore does impact the correlation and may 
be responsible for a dispersion in the direction orthogonal to the 
correlation, but according to the opacities, it cannot solely ex- 
plain the observed trend. 



5.3.3. Hypothesis for the depletion of submicron-sized grains 

Several solutions could be envisaged to explain the apparent de- 
pletion of submicron-sized grains at high altitudes above the 
disk midplane. The fragmentation time-scale of grains in disk 
atmospheres may be longer than the grain growth timescale, 
thereby quickly eliminating the smallest and most freshly pro- 
duced submicronic grains. The high velocity fluctuations ob- 
served in MHD simulations at high altitudes, do not, in prin- 
ciple, work in favor of this scenario. An alternative possibility 
may be that the submicron-sized grains are evacuated from the 
upper layers of the disks by stellar winds or radiation pressure. 
The impact of ra diation pressure on gra ins in young disks has 
been explored bv lTakeuchi & Liril (|2003|) and present several ad- 
vantages. First, it acts on short timescales as long as the grain 
size reaches a threshold, usually referred to as the blowout size 
limit, below which the radiation pressure force overcomes the 
gravitational and gas drag forces. Second, the blowout size limit 
falls in the right regime of grain sizes (close to a micrometer for 
silicates) for solar-type stars in the presence of some gas. 

It is interesting to apply the l atter scenario to the case o f 
more luminous stars. According to lKessler-Silacci et al.l ( l2007h . 
HAeBe stars tend to show a widest variety of 10 pm siUcate 
features, including the peaked features attributed to submicron- 
sized grains. Because of their higher temperature and luminos- 
ity, HAeBe stars are expected to have upper layers that are even 
more ionized than those of T Tauri stars, thereby enhancing the 
turbulence if it is supported by disk magnetic fields, as well as 
velocity fluctuations. This would therefore encourage the sys- 
tematic production of submicron dust fragments in disk atmo- 
spheres. On the other hand, silicate grains of a few micrometers 
and smaller, are in principle more prone to be eliminated by ra- 
diation pressure. Nevertheless, the latter reasoning assumes an 
optically thin medium at the wavelengths where grains absorb 
most, namely in the UV and visible. The optically thin region 
at mid-lR wavelengths is likely to be essentially opaque in the 
UV/visible, and this opacity will affect more strongly the most 
luminous stars, thereby possibly reducing the blowout size limit. 

Clearly, the depletion of submicron-sized grains in the disk 
atmospheres of most TTauri stars should be further investi- 
gated by models to examine the likelihood of the scenarios dis- 
cussed above (competition between coagulation and fragmenta- 
tion, stellar winds, radiation pressure, ...) 
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5.4. On radial mixing 

The coldest disk component probed by our IRS observations 
shows a rather high crystallinity fraction. According to Fig.|3] 
the 23 and 28 yum complexes are detected in more than 50% of 
the spectra. In a scenario where silicate crystals form by ther- 
mal annealing, the production of crystalline silicates can only 
occur on reasonably short time-scales near the star (T > 900 K), 
namely in regions much hotter than those emitting at /I > 20 fim. 
The observed cold crystals would thus have been radially trans- 
ported there from the inner disk zone, which according to two- 
dimension al models can occur around the midp lane of accre- 
tion d isks ( ICieslall2009l and references therein). iKeller & Gaill 
(120041) . for example, find that ~ 5 to 30% of the total mass that is 
transported inward at high disk altitudes is transported outward 
near the disk midplane. For a s olar mass star with an accretion 
rate of M = lO'^Mf^^/vr. lCieslal ([2009) estimates the crystalline 
fraction at 10 AU to be about ~ 40% by the end of his simula- 
tions. As the accretion rate becomes smaller, the crystalline frac- 
tion drops accordingly: ~ 10% and ~ 1% for M - lO'^M p/yr 
and 10 '^Mo/yr, respectively. However. IWatson et al.l ( |2009|) find 
only a weak correlation between the accretion rate and the crys- 
talline features arising at around lOfim, and no correlation with 
the 33 fim feature, leading the authors to suggest that some 
mechanisms may be erasing such trends. Alternatively, other 
crystallisa tion processes at colder temperatures (e.g. via nebu- 
lar shock, iDesch & Connollvll2002l) could also contribute to the 
presence of crystals at signifi cant distances from the central ob- 
ject (e.g. lKimura et al.ll2008h . 

It is noteworthy that the warmer disk zone seems signifi- 
cantly less crystalline than do the colder regions, although de- 
tailed compositional analysis of the IRS spectra is required 
before firmly concluding this must be so (Olofsson et al. 
2009b, in prep.) Interest ingly, composition al fits to some in- 
dividu al IRS spectra by iBouv et al.1 ( l2008h and ISargent et al] 
( |2009|) do show the need for larger amounts of silicate crys- 
tals in the outer regions compared to inner regions (e.g. 2MASS 
J04442713H-2512164and F04147+2822). If radial mixing leads 
to the large crystalline fraction in the cold disk component, it 
leaves the disks with an inhomogeneous chemical composition 
for what concerns the solid phase. Some re-amorphisation pro- 
cessing of the crystalline grains, by X-ray emission or cosmics 
rays for example, could also contribute to the apparent lack of 
crystals in the inner regions of disks. 

This result of a cold component being more crystalline com- 
pared to the inner warm regions see ms, at first glance , to be 
in contradiction with the results by Ivan Boekel et alj (12004 
who observed three HAeBe stars with the mid-IR instrument 
MIDI for the VLTI interferometer Using different baselines, 
they found that the inner regions (1-2 AU) of the disks are more 
crystalline compared to the outer regions (2-20 AU). Beside the 
fact that their sample is limited to 3 intermediate mass objects 
while we analyse a statistically significant TTs sample, direct 
comparison of their observations with ours is not straightfor- 
ward. They do have spatial information, but only for the lOyum 
spectral region, while we cover a larger spectral range. We, on 
the other hand, are observing the entire disk with Spitzer/IRS, 
with no direct spatial information. An ideal and challenging way 
to confront such kinds of interferometric data with IRS would be 
to repeat similar observations for TTs and even better, at wave- 
lengths larger than 20yum. This will ultimately require space- 
born interferometers because of the poor atmosphere transmis- 
sion in the mid-IR. With such observations, however, we would 



then be able to constrain the crystalline fraction in the inner and 
outer regions of the disk more directly. 

5.5. On tlie similarities with Solar System objects 

Studies of Solar Sytem asteroids and comets show strong sim- 
ilarities with our results. First of all, laboratory studies of as- 
teroidal and cometary solids have show that they contain very 
Uttle ISM-like materials. The SlPAVi ld 2 samples from the 
Stardust mission lMcKeegan et al.l ( |2006|) . for example, show that 
this comet contains high-temperature silicates and oxide miner- 
als. These authors conclude that such materials could not have 
been formed via annealing of presolar amorphous silicates in 
the Kuiper belt, and also that they could not have been formed 
from a single isotopic reservoir Overall, this suggests that comet 
8 IP/Wild 2 sampled different regions of the inner solar pro- 
toplanetary disk during its formation. One possible transporta- 
tion mechanism proposed is the combination of winds associ- 
ated with bipolar outflows. This hypothesis is also s upported 
by laboratory measurements from lToppani et all ( l2006h on con- 
densation under high temperature, low-pressure conditions. Mg- 
rich silicates crystals can rapidly (~ Ih at 4 x 10""' bar) be 
produced from the gas phase through condensation, confirm- 
ing the fact that crystalline grains can even be produced in 
bipolar outflows from evolved s tars. In still fu r ther w ork on 
the 8 IP/Wild 2 Stardust mission, IZolenskv et all ( |2006|) . found 
large abundances of crystalline materials. Olivine is present in 
most of the studied particles, with grain sizes ranging between 
sub micron-sized to over lOyum, and it is concluded that some 
materials in the comet have seen temper atures possibly higher 
than 2000 K. Similar to lMcKeegan et al.l (120061) . large-scale ra- 
dial transportation mechanisms inside the disk are required to 
match the derived mineralogy of the olivine samples. Similar 
results were found studying other solar system bodies, for ex- 
emple, comet Hale-Bopp. Spectroscopic studies of this object 
d Wooden et alJll999llWooden et al.ll2000l) revealed the presence 
of crystalline olivine, crystalline ortho-pyroxene and crystalline 
clino-pyroxene, in significant quantities. 

All these results echo the crystallinity paradox we find in our 
Spitzer observations, in the sense that the disk dust content is far 
from being homogeneous, and even far from the simple picture 
with crystalline grains close to the central object and amorphous 
content in the outer regions. 

6. Summary and conclusion 

We have conducted in this paper a comprehensive statistical 
study of crystalline silicates in proto-planetary (Class II) disks 
which was made possible thanks to the high sensitivity of the 
spectroscopic instrument IRS (5-35 fim) onboard the Spitzer 
Space Telescope. This unprecedented work is among the first 
statistical studies of crystalline silicates emission features over 
a large sample of young solar-analog stars (96 objects). We find 
the following results; 

1 . Crystallinity is not a marginal phenomenon in Class II disks 
around young solar analogs, as more than 3/4 of the ob- 
jects show at least one crystalline feature in the IRS spec- 
tra. Crystalline features emitting at wavelengths larger than 
20 fim are in addition widely present: both the C23 and C28 
complexes, produced by enstatite plus forsterite features, are 
present at a frequency larger than 50%. 

2. The crystalline silicates revealed by the IRS spectra are es- 
sentially Mg-rich silicates (forsterite, enstatite), with a small 
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Fig. 14. Schematic view of the 
TTauri disk regions probed by 
the Spitzer IRS spectroscopic 
observations. This illustrates 
some of the dynamical processes 
that are discussed in Sec.|5] and 
summarized in Sec.|6l 
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contribution of CaMg-rich silicates (diopside). We find no 
evidence of Fe-rich silicates (marginal evidence for fayalite) 
nor troilite (FeS) in our spectra. 

3. We find that the amorphous 10 //m feature is not correlated 
with the C23, C28 and 33.6/im crystalline features, neither 
on the basis of the frequencies at which they are detected 
nor on the fluxes they emit. Crystalline features emitted at 
wavelengths larger than 20 /^m are, on the other hand, corre- 
lated to each other (C23, C28 and 33.6 //m), based either on 
detection frequencies and band fluxes. 

4. From previous results, we conclude that within the spec- 
tral range of the Spitzer/IRS instrument, we are essentially 
probing two independent dust populations (see the illustra- 
tion in Fig.[T4b: a warm component, which is likely located 
close to the central object (< 1 AU), emitting features in 
the 10 /im region, and a second component, located at larger 
radii (< 10 AU) or deeper inside the disk; the latter cold com- 
ponent being responsible for the emission of both C23 and 
C28 complexes plus the 33.6yL(m forsterite feature. 

5. For a large majority of objects, the upper layers of the disks 
are essentially populated by micron-sized grains, suggest- 
ing that submicron grains are largely depleted in disks at- 
mospheres. For the amorphous lO/zm feature, we find the 
known correlation between the shape and the strength of 
the feature used as a grain size proxy and showing that the 
grains are /im-sized or larger. We adapted this relationship to 
the C23 complex and find similar results for the crystalline 
grains producing this complex. We also noted that the grain 
size proxies for the 10 fim and C23 features are uncorrected. 

6. We find a correlation between the shape of the SED 
(Fio/Fi^) and the grain size proxy for the 10//m feature. 
We see rather small, warm amorphous silicates for flared 
disks and grains becoming larger as disks are flattening. This 
trend is almost reversed for the cold, C23 crystalline grains 
as small crystals are only found for flattened SEDs. 

7. We identify an apparent crystallinity paradox since the cold 
crystalline features {A > 20 jum) are much more frequently 
detected than the warm ones (A ~ lOyum). We show that 
crystalline features can be hidden by amorphous features, 
especially near A ^ lOjjm, but this contrast effect is not 
sufficient to explain the difference regarding the apparition 
frequencies of the C23 complex and the 11.3yL(m forsterite 



feature, suggesting that the warm grain population is, on av- 
erage, intrinsically more amorphous than the cold grain pop- 
ulation. 

We argue that our Spitzer spectroscopic observations of disks 
around young stars reveal ongoing dynamical processes in planet 
forming regions of disks that future theoretical disk evolution 
and planet formation models will have to explain: 

- some mechanisms must supply the upper layers of the disks 
with micron-sized grains that should otherwise grow and 
settle rapidly. Vertical turbulent diffusion, accompanied by 
grain-grain fragmentation, would be natural candidates, 

- some mechanisms must act efficiently to remove the 
submicron-size grains from the regions which are optically 
thin at IRS wavelengths. Stellar winds and radiation pressure 
have been suggested as possible processes, 

- and, some mechanisms are responsible for rendering the 
colder disks regions much more crystalline than the warmer 
zones (crystallinity paradox). If this result has to do with 
dust radial mixing, it teaches us that the mixing does not re- 
sult in a disk with a homogenous grain composition, and that 
it must transport crystals efficiently into the "comet-forming 
regions". 

Acknowledgements. The authors warmly thank Bram Acke for his help on the 
IDL routine to extract the crystalline features characteristics, Jerome Aleon and 
Matthieu Gounelle for sharing with us their precious knowledge of small solid 
Solar System bodies, and Joel Green for the rich discussions we had. We also 
thank the anonymous referee for the very useful comments that helped improv- 
ing this study. It is also a pleasure to acknowledge fruitful discussions with the 
Grenoble FOST team members, as well as the participants to the ANR (Agence 
Nationale de la Recherche of France) project ANR-07-BLAN-0221. This re- 
search is partly supported by the Programme National de Physique Stellaire 
(PNPS). 

References 

Acke, B. & van den Ancker, M. E. 2004, A&A, 426, 151 

Alcala, J. M., Krautter, J., Schmitt, J. H. M. M., et al. 1995, A&AS, 114, 109 

Andrews, S. M. & Williams, J. R 2007, ApJ, 659, 705 

Bouwman, J., Henning, T., Hillenbrand, L. A., et al. 2008, ApJ, 683, 479 

Bouwman, J., Meeus, G., de Koter, A., et al. 2001, A&A, 375, 950 

Bouy, H., Huelamo, N., Pinte, C., et al. 2008, A&A, 486, 877 

Brown, J. M., Blake, G. A., Dullemond, C. R, et al. 2007, ApJ, 664, L107 

Chen, H., Grenfell, T. G., Myers, R C., & Hughes, J. D. 1997, ApJ, 478, 295 



Olofsson et al.: Crystalline silicates in TTauri disks 



19 



Chihara, H., Koike, C, Tsuchiyama, A., Tachibana, S., & Sakamoto, D. 2002, 

A&A, 391, 267 
Ciesla, F. J. 2009, Icarus, 200, 655 

Crovisier, J., Leech, K., Bockelee-Morvan, D., et al. 1997, Science, 275, 1904 
D'Alessio, P., Calvet, N., Hartmann, L., Franco-Hernandez, R., & Servm, H. 

2006, ApJ, 638,314 

Davoisne, C, Djouadi, Z., Leroux, H., et al. 2006, A&A, 448, LI 
Desch, S. J. & Connolly, Jr, H. C. 2002, Meteoritics and Planetaiy Science, 37, 
183 

Dominik, C, Dullemond, C. P., Waters, L. B. F M., & Walch, S. 2003, A&A, 
398, 607 

Dorschner, J., Begemann, B., Henning, T., Jaeger, C, & Mutschke, H. 1995, 

A&A, 300, 503 
Dullemond, C. P & Dominik, C. 2005, A&A, 434, 971 
Dullemond, C. P & Dominik, C. 2008, A&A, 487, 205 
Evans, N. J., Dunham, M. M., J0rgensen, J. K., et al. 2009, ApJS, 181, 321 
Evans, H, N. J., Allen, L. E., Blake, G. A., et al. 2003, PASP 1 15, 965 
Fabian, D., Henning, T., Jager, C, et al. 2001, A&A, 378, 228 
Fromang, S. & Nelson, R. P 2009, A&A, 496. 597 
Furlan, E., Hartmann, L., Calvet, N., et al. 2006, ApJS, 165, 568 
Gail, H.-P 1998, A&A, 332, 1099 
Gail, H.-P 2004, A&A, 413, 571 

Geers, V. C. 2007, PhD thesis, Leiden Observatory, Leiden University, P.O. Box 

9513, 2300 RA Leiden, The Netherlands 
Geers, V. C, Augereau, J.-C, Pontoppidan, K. M., et al. 2006, A&A, 459, 545 
Hallenbeck, S. L., Nuth, III, J. A., & Nelson, R. N. 2000, ApJ, 535, 247 
Hai-ker, D. E. & Desch, S. J. 2002, ApJ, 565, L109 

Harker, D. E., Woodward, C. E., Wooden, D. H., Fisher, R. S., & Trujillo, C. A. 

2007, Icai-us, 190, 432 

Herbig, G. H. & Bell, K. R. 1988, Third catalog of emission-line stars of the 

Orion population., ed. G. H. Herbig & K. R. Bell 
Honda, M., Kataza, H., Okamoto, Y. K., et al. 2003, ApJ, 585, L59 
Honda, M., Kataza, H., Okamoto, Y. K., et al. 2006, ApJ, 646, 1024 
Hughes, J. & Hartigan, P 1992, AJ, 104, 680 
Jaeger, C, Molster, E J., Dorschner, J., et al. 1998, A&A, 339, 904 
Keller, C. & Gail, H.-R 2004, A&A, 415, 1177 
Keller, L. P, Hony, S., Bradley, J. P, et al. 2002, Nature, 417, 148 
Kessler-Silacci, J., Augereau, J.-C, Dullemond, C. P, et al. 2006, ApJ, 639, 275 
Kessler-Silacci, J. E., Dullemond, C. P, Augereau, J.-C, et al. 2007, ApJ, 659, 

680 

Kimura, Y., Miyazaki, Y, Kumamoto, A., Saito, M., & Kaito, C. 2008, ApJ, 680, 
L89 

Koike, C, Chihara, H., Tsuchiyama, A., et al. 2003, A&A, 399, 1101 
Koike, C, Tsuchiyama, A., Shibai, H., et al. 2000, A&A, 363, 1115 
Ki-autter, J., Wichmann, R., Schmitt, J. H. M. M., et al. 1997, A&AS, 123, 329 
Lahuis, F, Kessler-Silacci, J. E., Evans, N. J., I., et al. 2006, c2d Spectroscopy 

Explanatory Supplement, Tech. rep., Pasadena: Spitzer Science Center 
Lahuis, R, van Dishoeck, E. E, Blake, G. A., et al. 2007, ApJ, 665, 492 
Laibe, G., Gonzalez, J.-F, Fouchet, L., & Maddison, S. T. 2008, A&A, 487, 265 
Lawson, W. A., Feigelson, E. D., & Huenemoerder, D. P 1996, MNRAS, 280, 

1071 

McKeegan, K. D., Aleon, J., Bradley, J., et al. 2006, Science, 314, 1724 
Merin, B., Augereau, J.-C, van Dishoeck, E. R, et al. 2007, ApJ, 661, 361 
Min, M., Hovenier, J. W., & de Koter, A. 2005, A&A, 432, 909 
Min, M., Hovenier, J. W., Waters, L. B. R M., & de Koter, A. 2008, A&A, 489, 
135 

Molster, F J., Waters, L. B. F M., & Tielens, A. G. G. M. 2002, A&A, 382, 222 

Natta, A., Meyer, M. R., & Beckwith, S. V. W. 2000, ApJ, 534, 838 

Natta, A., Testi, L., Calvet, N., et al. 2007, in Protostars and Planets V, ed. 

B. Reipurth, D. Jewitt, & K. Keil, 767-781 
Nuth, J. A. & Johnson, N. M. 2006, Icanjs, 180, 243 
Oliveira, I., Merin, B., Pontoppidan, K. M., et al. 2009, ApJ, 691, 672 
Padgett, D. L.. Cieza, L., Stapelfeldt, K. R., et al. 2006, ApJ, 645, 1283 
Pinte, C, Padgett, D. L., Menai'd, F, et al. 2008, A&A, 489, 633 
Press, W. H., Teukolsky, S. A., Vetterling, W. T, & Flannery, B. P 1992, 

Numerical recipes in C. The art of scientific computing, ed. W. H. Press, 

S. A. Teukolsky, W. T. Vetterling, & B. P. Flannery 
Przygodda, R, van Boekel, R.. Abraham, P, et al. 2003, A&A, 412, L43 
Salyk, C, Pontoppidan, K. M., Blake, G. A., et al. 2008, ApJ, 676, L49 
Sargent, B. A., Fon-est, W. J., Tayrien, C, et al. 2009, ApJS, 182, 477 
Savitzky, A. & Golay, M. J. E. 1964, Analytical Chemistry, 36, 1627 
Servoin, J. & Piriou, B. 1973. Phys. Stat. Sol. (b), 55, 677 
Sitko, M. L., Lynch, D. K., & Russell, R. W. 2000, AJ, 120, 2609 
Sloan, G. C, Kraemer, K. E., Goebel, J. H., & Price, S. D. 2003, ApJ, 594, 483 
Takeuchi, T. & Lin, D. N. C. 2003, ApJ, 593, 524 

Toppani, A., Libourel, G., Robert, F., & Ghanbaja, J. 2006, 

Geochim. Cosmochim. Acta, 70, 5035 
Valenti, J. A., Fallon, A. A., & Johns-Krull, C. M. 2003, ApJS, 147, 305 



van Boekel, R., Min, M., Leinert, C, et al. 2004, Nature, 432, 479 
van Boekel, R., Min, M., Waters, L. B. F M., et al. 2005, A&A, 437, 189 
van Boekel, R., Waters, L. B. R M., Dominik, C, et al. 2003, A&A, 400, L21 
Watson, D. M.. Leisenring, J. M., Purlan, E., et al. 2009, ApJS, 180, 84 
Wichmann, R., Krautter, J., Schmitt, J. H. M. M., et al. 1996, A&A, 312, 439 
Wooden, D., Desch, S., Harker, D., Gail, H.-P, & Keller, L. 2007, in Protostai's 

and Planets V, ed. B. Reipurth, D. Jewitt. & K. Keil, 815-833 
Wooden, D. H., Butner, H. M., Haiter, D. E., & Woodwai'd, C. E. 2000, Icarus, 

143, 126 

Wooden, D. H., Harker, D. E., Woodward, C. E., et al. 1999, ApJ, 517, 1034 
Zolensky, M. E., Zega, T. J., Yano, H., et al. 2006, Science, 314, 1735 



Table 4. Characteristics of thie observations 



Starname 


RA (J2000) 


Dec (J2000) 


Class 


SpT 


AOR Key 


Obs. date 


Modules 


Perseus 


RN015 


03h27m47.68s 


30dl2m04.3s 


TTs (II) 


- 


0005633280 


2004-08-30 


SL, SH, LH 


LkHo' 325 


03h28m52.22s 


30d45m05.5s 


CTTs (II) 




0015737600 


2005-09-05 


SL, LL 


LkHa'270 


03h29ml7.66s 


31d22m45.1s 


TTs (II) 


K2.5-K7 


0005634048 


2005-02-08 


SL, SH, LH 


LkHQ;271 


03h29m21.87s 


31dl5m36.3s 


TTs (II) 


K3-K5 


0011827968 


2005-02-10 


SL, SH, LH 


SSTc2dJ033036.0+303024 


03h30m35.94s 


30d30m24.4s 


II 


- 


0015737600 


2005-09-05 


SL, LL 


SSTc2dJ033037.0+303128 


03h30m36.98s 


30d31m27.7s 


II 


- 


0015737600 


2005-09-05 


SL, LL 


LkHa 326 


03h30m44.05s 


30d32m46.7s 


TTs (II) 


G-MO 


0005634304 


2005-02-10 


SL, SH, LH 


SSTc2dJ033052.5+305418 


03h30m52.58s 


30d54ml8.1s 


11 


- 


0015737088 


2005-09-05 


SL, LL 


SSTc2dJ033241. 7+31 1046 


03h32m41.76s 


31dl0m46.6s 


II 




r\r\ 1 CTOTO A A 

0015737344 


"^rvrxc r\c\ f\c 

2005-09-05 


SL, LL 


LkHa 327 


03h33m30.41s 


31dl0m50.4s 


TTs (II) 


A9-K2 


0005634560 


2004-02-03 


SL, SH, LH 


SSTc2dJ033341. 3+31 1341 


03h33m41.30s 


31dl3m41.6s 


II 




0015736832 


2005-09-05 


SL, LL 


BD+31634 


03h41m39.20s 


31d36ml0.3s 


HAeBe (II) 


A8V 


0015736320 


2005-09-05 


SL, LL 


SSTc2dJ034219.3+314327 


03h42ml9.31s 


31d43m26.6s 


n 




0015736064 


2005-09-05 


SL, LL 


IRAS 03406+3133 


03h43m44.54s 


31d43m09.1s 






0015735296 


2005-09-05 


SL,LL 


LkHa 330 


03h45m48.29s 


32d24mll.8s 


CTTs (II) 


G3 


0005634816 


2004-09-02 


SL, SH, LH 


IRAS 03446+3254 


03h47m47.12s 


33d04m03.4s 


TTs (II) 




0005635072 


2004-09-29 


SL, SH, LH 


Taurus 


LkCaS 


04h24m57.08s 


27dllm56.5s 


CTTs (II) 


MO 


0009832960 


2005-02-08 


SH, LH 


IQTau 


04h29m5 1.56s 


26d06m45.0s 


WTTs (II) 


M0-M0.5 


0009832704 


2005-02-09 


SH,LH 


FXTau 


04h30m29.62s 


24d26m45.1s 


C+WTTs (II) 


M1-M4 


0009832448 


2005-02-10 


SH,LH 


V710Tau 


04h31m57.79s 


18d21m36.3s 


C+WTTs (II) 


M0.5-M3 


0005636608 


2004-09-29 


SH,LH 


RXJ0432.8+1735 


04h32m53.23s 


17d35m33.7s 


TTs (II) 


M2* 


0015917824 


2005-09-09 


SL,LL 


DNTau 


04h35m27.37s 


24dl4m58.8s 


CTTs (II) 


MO 


0009831936 


2005-02-10 


SH,LH 


CoKu Tau 3 


04h35m40.94s 


24dllm08.6s 


CTTs (II) 


Ml 


0009831936 


2005-02-10 


SH,LH 


CoKuTau4 


04h41ml6.79s 


28d40m00.5s 


CTTs (II) 


M1.5 


0005637888 


2004-09-02 


SH,LH 


RRTau 


05h39m30.52s 


26d22m27.0s 


HAeBe (II) 


B8-A5 


0005638400 


2004-09-28 


SL, SH, LH 


Chamaleon 


SXCha 


10h55m59.74s 


-77d24m39.9s 


TTs (II) 


M0.5 


0005639424 


2004-08-31 


SH,LH 


SYCha 


10h56m30.47s 


-77dllm39.4s 


TTs (II) 


MO 


0005639424 


2004-08-31 


SH,LH 


TWCha 


10h59m0 1.10s 


-77d22m40.80s 


TTs (II) 


KO-MO 


0005639680 


2004-09-01 


SH,LH 


B35 


Ilh07m2 1.49s 


-77d22mll.8s 


CTTs (II) 


M2 


0005639680 


2004-09-01 


SH,LH 


VWCha 


Ilh08m01.50s 


-77d42m28.70s 


CTTs (II) 


K2-K7 


0005639680 


2004-09-01 


SH,LH 


Hn9 


Ilh09ml8.16s 


-76d30m29.2s 




K4-M1'' 


0009831168 


2005-02-10 


SL, LL 


VZCha 


Ilh09m23.80s 


-76d23m20.7s 


TTs (II) 


K6-K7 


0005640448 


2004-09-02 


SH,LH 


WXCha 


Ilh09m58.75s 


-77d37m08.9s 


TTs (II) 


K7-M0 


0005640192 


2004-09-01 


SH,LH 


ISO-Cha237 


llhlOml 1.44s 


-76d35m29.2s 


TTs (II) 


MO 


0005640448 


2004-09-02 


SL,LL 


C7-11 


Ilhl0m38.01s 


-77d32m39.9s 


TTs (II) 


KB 


0005640192 


2004-09-01 


SH,LH 


HM 27 


Ilhl0m49.62s 


-77dl7m51.7s 


TTs (II) 


K7 


0005640192 


2004-09-01 


SH,LH 


XXCha 


llhllm39.67s 


-76d20ml5.1s 


TTs (II) 


M1-M2 


0005640448 


2004-09-02 


SH,LH 


HD 98922 


Ilh22m3 1.67s 


-53d22mll.4s 


HAeBe (II) 


B9 


0005640704 


2004-01-04 


SH,LH 


HD 101412 


Ilh39m44.46s 


-60dl0m27.7s 


HAeBe (II) 


B9.5 


0005640960 


2005-02-10 


SL, SH, LH 


TCha 


Ilh57ml3.53s 


-79d21m31.5s 


TTs (II) 


G2-K0 


0005641216 


2004-07-18 


SH,LH 


IRAS 12535-7623 


12h57ml 1.78s 


-76d40mll.5s 


TTs (II) 


MO 


0011827456 


2004-08-31 


SH,LH 



Table 4. continued. 





RA ri2000~) 


Dec ("12000^ 


Class 


SpT 


AOR Kev 






ISO-Chall 13 


12h58m06.70s 


-77d09m09.5s 


BD (II) 




0015918592 


2005-09-13 


SL, LL 


RXJ1301.0-7654 


13h00m53.23s 


-76d54ml5.2s 




Kl'' 


0009830656 


2005-03-12 


SL,LL 


Sz50 


13h00m55.37s 


-77dl0m22.2s 


TTs (II) 


K7-M3 


0011827456 


2004-08-31 


SH,LH 


ISO-ChaII54 


13h00m59.20s 


-77dl4m02.7s 


TTs (II) 


K3 


0015735040 


2005-08-13 


SL, SH, LL 


Sz52 


13h04m24.90s 


-77d52m30.3s 






0018215168 


2006-07-03 


SL, LL 


Sz62 


13h09m50.66s 


-77d57m24.0s 




wr 


0009830400 


2005-03-12 


SL, LL 


Lupus 


HTLup 


15h45ml2.87s 


-34dl7m30.6s 


TTs (II) 


K2 


0005643264 


2004-08-30 


SL, SH, LH 


GWLup 


15h46m44.68s 


-34d30m35.4s 


TTs (II) 


M2-M4 


0005643520 


2004-08-30 


SL, SH, LH 


HM Lup 


1 CU/1*7™CA iCO^ 

Ijn4/mj0.o3s 


-3jd/am3j.4s 


lis (11) 


N /T /I 

M4 


0005643 / /6 


/004-0a-30 


CT T T 

SiL, LL 


Sz73 


15h47m56.98s 


-35dl4m35.1s 


TTs (II) 


K2-M 


0005644032 


2004-08-30 


SL, SH, LH 


GQLup 


15h49ml2.10s 


-35d39m05.0s 


TTs (II) 


K7-M0 


0005644032 


2004-08-30 


SL, SH, LH 


Sz76 


15h49m30.74s 


-35d49m51.4s 




Ml-'' 


0015916288 


2005-09-09 


SL, LL 


IMLup 


15h56m09.20s 


-37d56m06.4s 


TTs (II) 


MO 


0005644800 


2004-08-30 


SL, SH, LH 


RULup 


15h56m42.31s 


-37d49ml5.5s 


CTTs (II) 


K3-M0 


0005644800 


2004-08-30 


SL, SH, LH 


Sz84 


15h58m02.50s 


-37d36m02.8s 




M5.5« 


0005644288 


2004-03-25 


SL, LL 


RYLup 


15h59m28.39s 


-40d21m51.2s 


TTs (II) 


K0-K4 


0005644544 


2004-08-30 


SL, SH, LH 


EX Lup 


16h03m05.50s 


-40dl8m24.9s 


TTs (II) 


MO 


0005645056 


2004-08-30 


SL, SH, LH 


RXJ1603.2-3239 


16h03ml 1.81s 


-32d39m20.2s 


TTs (II) 


K7« 


0015917312 


2005-09-09 


SL,LL 


Sz96 


16h08ml2.62s 


-39d08m33.4s 




Ml.5^ 


0016755200 


2006-03-15 


SL, LL 


Szl02 


16h08m29.70s 


-39d03mll.2s 


TTs (II) 


MO 


0009407488 


2004-03-25 


SL, SH, LH 


SSTc2dJ161 159.8-382338 


16hllm59.81s 


-38d23m37.5s 


II 




0015737856 


2005-08-14 


SL, LL 


RXJ1615. 3-3255 


16hl5m20.23s 


-32d55m05.1s 


TTs (II) 


K5'' 


0015916800 


2005-09-09 


SL, LL 


Ophiuchus 


AS 205 


16hllm31.35s 


-18d38m26.1s 


TTs (II) 


K5 


0005646080 


2004-08-28 


SL, SH, LH 


Harol-1 


16h21m34.69s 


-26dl2m27.0s 


CTTs (II) 


K5-K7 


0009833472 


2005-03-12 


SL, SH, LH 


SSTc2dJ162148.5-234027 


16h21m48.49s 


-23d40m27.4s 


II 




0015920897 


2006-03-15 


SL,LL 


SSTc2dJ 1 6222 1 .0-230403 


16h22m21.01s 


-23d04m02.6s 


I 




0015920897 


2006-03-15 


SL, LL 


SSTc2dJ162245.4-243124 


16h22m45.40s 


-24d31m23.9s 


II 




0015920641 


2006-03-15 


SL, LL 


SSTc2dJ162332.8-225847 


16h23m32.85s 


-22d58m46.9s 


II 




0015920641 


2006-03-15 


SL,LL 


Harol-4 


16h25ml0.51s 


-23dl9ml4.5s 


TTs (II) 


K4-K6 


0009833216 


2005-03-12 


SH,LH 


VSSGl 


16h26ml8.86s 


-24d28ml9.7s 


TTs (II) 




0005647616 


2004-08-28 


SH,LH 


DoAr24E 


16h26m23.38s 


-24d21m00.1s 


TTs (II) 


KO-Kl 


0005647616 


2004-08-28 


SH,LH 


Do At 25 


16h26m23.68s 


-24d43ml4.00s 


TTs (II) 


K5'' 


0012663808 


2005-08-12 


SH,LH 


GY23 


16h26m24.06s 


-24d24m48.1s 


TTs (II) 


K5-M2^' 


0005647616 


2004-08-28 


SH,LH 


SR21 


16h27ml0.28s 


-24dl9ml2.5s 


TTs (II) 


F4-G2.5 


0005647616 


2004-08-28 


SH,LH 


SSTc2dJ162715. 1-245139 


16h27ml5.14s 


-24d51m38.9s 


II 




0016754688 


2006-03-15 


SL,LL 


SR9 


16h27m40.27s 


-24d22m04.0s 


TTs (II) 


K5-M2 


0012027392 


2004-09-02 


SH, LH 


SSTc2dJ162816.7-240514 


16h28ml6.73s 


-24d05ml4.3s 


II 




0016754944 


2006-03-15 


SL, LL 


V8530ph 


16h28m45.28s 


-24d28ml9.0s 


TTs (II) 


M1.5 


0012408576 


2005-03-12 


SH,LH 


ROX42C 


16h31ml5.75s 


-24d34m02.1s 


TTs (II) 


K4-K6 


0006369792 


2005-03-12 


SH,LH 


ROX43A 


16h31m20.13s 


-24d30m05.1s 


TTs (II) 


GO 


0015914496 


2005-09-09 


SH,LH 


IRS60 


16h31m30.89s 


-24d24m39.7s 






0006370048 


2005-09-09 


SH,LH 


Harol-16 


16h31m33.47s 


-24d27m37.1s 


TTs (II) 


K2-K3 


0012664064 


2005-09-09 


SH, LH 



Table 4. continued. 



Starname 


RA (J2000) 


Dec (J2000) 


Class 


SpT 


AOR Key 


Obs. date 


Modules 


Harol-17 


16h32m21.94s 


-24d42ml4.7s 


TTs (II) 


M2.5 


0011827712 


2004-08-29 


SL, SH, LH 


RNO90 


16h34m09.20s 


-15d48ml6.80s 


TTs (II) 


G5 


0005650432 


2004-08-28 


SL, SH, LH 


WaOph6 


16h48m45.63s 


-14dl6m35.96s 


TTs (II) 


K7 


0005650688 


2006-03-15 


SL, SH, LH 


V11210ph 


16h49ml5.31s 


-14d22m08.6s 


CTTs (II) 


K5 


0005650688 


2006-03-15 


SL, SH, LH 


HD 163296 


17h56m21.29s 


-21d57m21.9s 


HAeBe (II) 


A0-A2 


0005650944 


2004-08-28 


SH, LH 


Serpens 


V V ocr 


i onzoni^ / .oos 


uauonijy.oS 


ri/\eDe \\\-} 


U 1 A 






oij, orl, i^rl 


SSTc2d J 182850.2+00950 


18h28m50.21s 


0d09m49.6s 


II 


M2*' 


0013461505 


2006-04-21 


SL, LL 


SSTc2d J 1 82900.9+0293 1 


18h29m00.90s 


0d29m31.5s 


n 




00132101 12 


2005-04-17 


SL, SH, LL 


1 C/u J io/yUi .z+U/Vj J 


i on/ymu i ./us 


uci/ymj j./s 


ii 


MU.J 


UU ij40ijUj 


/UU0-U4-/1 


CT T T 


SSTc2dJ182901.8+02954 


18h29m01.80s 


0d29m54.3s 


n 


K7* 


0013210112 


2005-04-17 


SL, SH, LL 


SSTc2dJ182909.8+03446 


18h29m09.80s 


0d34m45.8s 




MS 


0013210624 


2005-04-14 


SL, SH, LL 


SSTc2d J 182928.2+02257 


18h29m28.24s 


-0d22m57.4s 


n 




0013210368 


2005-04-14 


SL, SH, LL 


EC69 


18h29m54.35s 


Idl5m01.8s 


II 




0009407232 


2004-03-27 


SL, LL 


EC82 


18h29m56.89s 


Idl4m46.5s 


TTs (II) 


MO 


0009407232 


2004-03-27 


SL, SH, LH 


EC90 


18h29m57.75s 


Idl4m05.9s 


TTs (II) 




0009828352 


2004-09-01 


SL, SH, LH 


EC92 


18h29m57.88s 


Idl2m51.6s 


TTs (II) 


K7-M2 


0009407232 


2004-03-27 


SL, SH, LH 


CK4 


18h29m58.21s 


Idl5m21.7s 


TTs (II) 


K3 


0009407232 


2004-03-27 


SL, SH, LH 


LkHo' 348 


18h34ml2.65s 


-0d26m21.7s 


TTs (II) 


F6 


0009831424 


2006-04-16 


SL,SH, LH 


BFOri 


05h37ml3.26s 


-6d35m00.6s 


HAeBe (II) 


A5-F6 


0005638144 


2004-10-03 


SL, SH, LH 


IRAS 08267-3336 


08h28m40.70s 


-33d46m22.3s 


TTs (II) 


K2-K3 


0005639168 


2004-11-11 


SL, SH, LH 


HD 135344 


15hl5m48.44s 


-37d09ml6.0s 


HAeBe (II) 


A0-F4 


0005657088 


2004-08-08 


SH,LH 



Where reference for spectral type is not indic ate d, it comes from Lahuis et al. ( 2007h. the thir d p aper in the series of four c2d ar ticles on IRS obse rvat ions of TTauri s tars, a ) 
iHerbig & Belli (Il988l) b Wichmann et al.l (Il996[). c'llTawson et alJ (IT996h. d')lAlcala et al.l (Il995h. e-jlHughes & HartiganI (Il992h. fi lChen et al.l(ll997h . g^l lKrautter et aP (Il997h . 
h) |Pad^ett et all ( l2006l) . D iAndrews & Williamsl ( |2007|) . i) lKessler-Silacci et all ( l2006l) . k) l01iveira et all ( l2009l) 



Table 5. Signal-to-noise ratios for detections of crystalline silicate features, 
amorphous silicate 10 /im feature and silica. 



Starname 


9.2 


Amorphous 10.0 


11.3 


12.5 


16.2 


21.6 


C23 


25.0 


C28 


33.6 


S 11.3/5 9.8 


Peak 


S 24 /S 23 


^ 23jUH7 

Peak 


F3o/Fn 


Perseus 


RN015 


- 


871 


- 


- 


- 


- 


- 


- 


- 


- 


0.92 


1.45 


- 


- 


1.13 


LkHa 325 


39 


599 


- 


- 


- 


- 


35 


- 


39 


40 


0.98 


1.51 


- 


- 


1.10 


LkHa270 


22 


109 


- 


- 


- 


- 


108 


- 


- 


- 


1.01 


1.11 


0.51 


9.83 


1.59 


LkHa271 


- 


106 


- 


13 


- 


- 


- 
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Fig. 15. Spitzer/IRS spectra of young stars 
estimated in Sec. 12. 41 Sources are sorted by 
and dark purple for LH. 



in the Perseus cloud, in units of Jy. Gray areas correspond to 3-cr uncertainties as 
increasing Right Ascension. Red color is for SL module, orange for SH, blue for LL 
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Fig. 16. Same as Fig.[T5]but for the Taurus star sample 
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Fig. 17. Same as Fig.[T5]but for the Chamaleon star sample. 
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Fig. 19. Same as Fig.[T5]but for the Lupus star sample 
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Fig. 20. Same as Fig.[T5]but for the Ophiuchus star sample. 

AS 205 Haro 1-1 



SSTc2d J16214B. 5-234027 






O.B - 




0.7 [ 








0.6 - 








0.5 r 




0.4 r 



10 15 20 25 30 35 
SSTcSd J162221. 0-230403 




10 15 20 25 30 35 
Haro 1-4 




10 15 20 25 30 35 
A [/Am] 

DoAr25 




10 15 20 25 30 35 
A [/^m] 

SSTc2d J162715. 1-245139 




10 15 20 25 30 35 
A [/im] 

V8530ph 




10 15 20 25 30 35 
A [/Lim] 




10 15 20 25 30 35 
A [/.im] 

SSTc2d J162245. 4-243124 



15 20 
A [/.tm] 

SSTc2d J162332.f 



25 30 35 
225847 




5 10 15 20 25 30 35 
A [fJ.m] 

VSSGl 



10 15 20 25 
A [/.im] 

DoAr 24E 




5 10 15 20 25 30 35 
A [/.im] 

GY23 




15 20 25 
A [/.im] 

SR 21 



35 





A [/Am] 
SR 9 



10 15 20 25 30 35 
A [/Am] 

SSTc2d J162816. 7-240514 





10 15 20 25 30 35 
A [/Am] 

R0X42C 



5 10 15 20 25 30 35 
A [/Am] 

R0X43A 





10 15 20 25 30 35 
A [/Am] 



5 10 15 20 25 30 35 
A [/Am] 



32 



Olofsson et al.: Crystalline silicates in TTauri disks 



Fig. 21. Spitzer/IRS spectra for the Ophiuchus sample: continued 
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Fig. 22. Same as Fig.[T5]but for the Serpens star sample. 
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Fig. 23. Same as Fig.[T5]but for the 3 isolated stars in our sample 
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